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Until recently there have been three standard iterative 
techniques for improving the accuracy of numerical solutions s 
Richardson extrapolation. Iterated Deferred Correction (IDC) 
of Fox and Poreyra, and the relatively new Iterated Defect 
Correction (IDcC). The method IDeC was originated by 
Zadunaisky as error estimation, and later formulated by 
Stettor for iterative apolication. The present th'^sis 
introduces now and generalized methods of iterative improve- 
ment and presents a variant of IDeC as well as various novel 
applications of IDeC. 

A historacal survey outlining the related developments 
is given in the introductory chapter of this thesis. 

Chapter 2 deals with now applications of IDeC based on 
quadrature methods for the solution of linear Fredholm and 
linear and nonlinear Volterra integral equations of the 
second kind as well as of linear eigenvalue problems. 



Chapters 3 to 6 present lOur now acculeratxon tcc'''n3aucs 
which haV'j been named as Successive Extrapolated IDoC 
( SEIDoC) , Iterated Defer icd Defect Correction (IDDcC), 
Successive Extrapolated IDDoC (SEIDDoC), and IDeC on 
E xtr apolo cion (IDeCE). These methods, which involve IDoC 
with tno use of i nt<. rpolnting polynomials, are dosciibod as 
applied to quadrature methods of Fredholm integral equations. 
The asymDtotic expansions for the methods are derived? 
implementation details arc provided? and illustrative nuinO" 
rical cxamolGs have boon worked out. Applications to linear 
and nonlinear Voi terra equations have also been considered. 

A variant oi iUeC railed derroq^ive Updated iDeC (SUIDoC) 
IS presented in Chapter 7. A conparotivr study of all those 
acceleration techniques is mace in Chaptoi h. Tho ninth 
chapter is devoted to applications of defect correction, with 
the use of splines, to different typos of partial differential 
equations. 

Tho class of acceleration techniques has been enriched 
remarkably by tho IDeC method, '//hen applied to Vol terra 
equations, it produces solutions of high accuracy oven for 
relatively large stopsizes. Fox Fredholm equations with 
smooth kernels and functions, tho easy-to-implement SEIDoC 
method extends tho order of convergence arbitrarily beyond 
the maximum order attainable by IDeC? each iteration of IDeC 
and IDDeC increases the order of convergence by n and n+p 



r e spectiveiv , where n is tne order of the basic discretization 
scheme and p is the jncr^^ase in oroer due to deferred corre- 
ct 3 on. Because of its goaranteed convergence, IDDoC is 
suponor to IDC also. The method of SEIDDcC extends the 
order of ccviv 'rgonce of IDDoC. 

numerical results confirm theoretical analysis and 
indicate the capability of the defect correction tochniquos 
in pr'^viding general software for the automatic computation 
of solutions to problems not only in integral equations but 
also in other operator equations. 

All the above methods have general applicability and v«/ill 
bo useful in various areas. The computations have been 
carried out using DEC 1090 at IIT, Kanpur and ICL 19043 at 


BHU, Varanasi. 



CHAPTER 1 


INT RODUCTIOR 


In tiiip thesis we discuss a numher of techniques for 
improving the accuracy of the numerical solution of operator 
equations of the form 

(1.1) RCy) = 0 

where the function y is the exact solution. 

In most of the discretization methods which are used 
for solving equation (l.l) the numerical solution at the 
discrete points converges to the true solution as th^ 
discretization parameter h t<^nds to zero. However the rate 
of convergence is often slow, and for satisfactory accuracy 
it becomes necessary to choose a very small h. Hence we 
need techniques to accelerate this rate of convergence or 
equivalently for procedures which estimate the error in 
the numerical solution with a view to improve the accuracy 
in the solution iteratively. During the past century a 
number of such techniques have been introduced in the 
liteirature. 

Since the beginning of this century, asymptotic 
expansions of the global error in powers of h have been 
used to construct methods of increasing accuracy. One of 
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the mrthods is that of RichardsorL[l9lO , 1927] called Toy 
him the 'deferred approach to the limit’. In this method, 
solutions of incr'=asing accuracy .are produced on the original 
grid by tho us^ of solutions on finer grids. 

Another well-known and ¥id''ly applicable acceleration 
t'^chniq.U'^ is th'^ 'difference correction' or 'deferred 
correction' method of I’oy[l947, 1949, 1950, 1957]. This 
method performs computations on the original grid of 
discr itisiation itself without refining the grid but uses 
high order approximations without increasing the complexity 
of tho linear or nonlinear system resulting from th"" dis- 
cretization scheme for (l.l). This method has later been 
analyzed, refined .and extended by Pereyr.a[1966, 1967a, 

1967b, 1968 , 1970 ], who named tho refinpd version of Fox’s 
ra<^thod 'It'^rated Deferred Correction' (IDC). Peri^yra 
used th.e r<^oults on the asymptotic expansions for the 
global error, such as in Gragg[l965] and Stetter [I 965 ] , and 
showed thnt for sufficiently small h each it'^nation of IDC 
improves the order of accuracy by p, where p is th"’ order 
of th(^ discretization scheme. The IDO method gives much 
better and faster r'^sults than most other methods of discre- 
tization in real pnactical problems (sec Fox [1977]). 

Recently an interesting and efficient approach to 
global ."'rror "estimation h.as been developed by tho Argen- 
tinian astronomer P.S. Zadunaisky. This new approach 



3 


discussed by Lambert [197 7j and the related method of 
Iterat'^d Defect Correction by Stetter [1974, 1978J end 
Frank[1976, 1977] stimulated the work presented in this 
thf^sis. Our primary goel was to opply Zadunaisky's de-vice 
iterat-vely in new ways and in new types of operator equa- 
tions other than ordinary differential equations. 

In the next section of this chapter historical out- 
lines of the development of Zaduneisky’s principle and the 
method of Iterated D-^foct Correction are givf^n. The r^-’la- 
t3onshi5)S bo tween this method and other methods of it-^rativo 
improvement, such as those introduced by Frank, Hortling and 
Ueberhub-^r end by Lindberg, are also pointed out, A summary 
of thc=' thesis IS given in Section 1*2. 

1 1 H istorical Survey 

Zadunaisky, when solving certain problems in Astro- 
nomy which invo]vi^d highly accurate observations found it 
impos.siblc or v-'ry difficult to apply the then kno-wn 
methods for the purpose of achieving tho desiri-^d accuracy 
in th<~' numerical solution of a system of ordinary differen- 
tial equntions ( Zadunaisky [1976] ) . So he was motivated by 
the desire to devise a good procedure for the estimation 
of th-' global errors propagated in the solution, and thereby 
to improve th? accuracy of that solution. 
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He successfully devf^lopod one such procedure and 
introduc<^d it first in a symposium of th-^ International 
Astronomic'^1 Union in 1964 (Zadunaisky [1966a] ) . Though he 
had also discussed this procedure in his subseluont 
papers ( Zadunaisky [l966hj 1970, 1972]) which deal 

with problems in astroromy, th'^ first time a group of 
numerical analysts came to know of his idea was at the 
Dunde'^ Conference on Numerical Analysis in 1973 when 
Zndunaisky presented a paper. In that paper he had included 
a mathematical description of his procedure, and impr'''ssiwe 
numerical results for problems in systems of differential 
•^q. nations, two-point boundary v-^lue problems, and appli- 
cation to lir.'^ar multist^p and predictor-corrector methods. 

St ittor [1974] subsequently found theoretical results 
juste-fying Zadumioky's numerical results. He formalized 
Zadunaisky 's idea of error estimation and conceived of its 
itentivo use. 

Th ’ pow ^r of ZrCdunaisky ’ F epproach is indicated in 
the math-^matical analysis of I'rank[l975] on the asymptotic 
b )haviour of thr approach when applied to Rung'^-Kutta 
methods of order p. Ho prov">d that under suitable smoothness 
conditions Zadunaisky’ s ostimatr' for global orror differs 
from th'^ actual global prror by OCh^^). 

Alf^ld[l975j gives a survjy of Zadunaisky ’s device 
applied to initial and boundary value problems of ordinary 
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dif f '^rsntial equations. The papers of Frank [1976, 1977] 
deal with iterativ-^ application of Zadunaisky ’ s concept 
of error estimation to two-point boundary value problems | 
he called this technique the method of ’Iterated Defect 
Correction' (iDeC), whereas Stetter [1974] had introduced 
the term Differential Correction in place of Defect Correc- 
tion. Moroov^r, Frank was the first to prove rigorous 
results concerning the asymptotic behaviour of IDeC applied 
to classical discretization methods for two-point boundary 
v-^lue problems. 

Frank and Ueberhuber[1975] proved that for an IDeC 
procedure based on an explicit or implicit Rungo-Kutta 
method of order p the sth iterate is of order s.p, the 
maximum attainable order of the procedure being limited by 
the degree of polynomials used in piecewise interpolation 
of the discrete solution. Besides, th-^ir IDeC applications, 
on the implicit Ruler scheme, have proved to be efficient 
in the solution of stiff systems of ordinary differential 
equations (Frank and Uebf=rhuber[l977] ) . 

Meanwhile, by adapting Zadunaisky 's idea Frank, 
Hertling and Ueberhuber[l976, 1977] introduced a n'^w tech- 
nique for estimating local discretization errors and a new 
variant to IDC methods based on that technique. Unlike 
the other IDC m'=^thods, these IDC methods do not depend on 
explicit knowledge of the asymptotic expansion for the 
local discretization error. Consequently, they have a 
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Inrgor dom^.an of ipplxcability coinpT.rod with other IDG 

impl-^Tii ntn,tiorL8 such as ttose of Perpyra[l973 ] and L-^ntini 
and P cr oyra [ 1975 ] , m which th'=> computation of deferred 
corr^c bicnn x& rabh'-^r difficult and is reliable only for 
£nDt’bj.'-' choices of the discretization parameter h (see 
Laniel end Martin[l977] ' . 

P'tirnulated by hadunaisky ' s approach discussed in 
St e : 'ce 3 ’[i 97 4 ] , lindb ;r^[l976] developed yet another tech- 
nicine Co-r ipepo__ improv 'meat. This technique involves 
" s'-qu'nce of perturbed problems but is quite differ 'nt 
Irou bhii of ZidrnaDSky. It if‘ ilso more general than the 
difference correction procedure of Pox[l957]. In Lindb erg's 
bechniquj Jocei error of.tame'-Qg err constructed avoiding 
d' *“'il^d loc J. error an'-Oysis. whereas in Per Lyra's method 
i 1 IS T'' '^jarary to use local Tror expansions. Lindborg 
po:'''’L.f' or+ ' ' 3 " advance e of Por'-yra’s version of IPG, 

’:h''t who 'tM'y 3 c nonlineer in some derivative of 
y„ che oxr)"'Lr LO}*s a,”o aerj, conbersome to find and ever more 
cumber 50m xo approxnia to ’ , Iincll rg also gives various 
exau.ples ^or th^ application of bis tachnique to ordinery 
and pertial diC Ceren'ciil equations end to integral aquations. 
7e shOvrs 'fiat for each iteration che order of accuracy oc 
the n-'w approximate solution is nc reused as in the cas"^ of 
Li)G. 

The notable difference between the IDeG method and 
Lindberg's technique is that in the former a sn o<oth global 
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ripproxiuntion to tho solution of tho operator equation 
F(y) =0 IS construct od an cxbBnsion of th^' discrete 
solution and th.m p'=-rturbation or def'^ct is comput'^d from 
that 'approximation, -whereas th^^ latter technique uses the 
discri^t'-' solution directly and the local properties of the 
p ••^rturhe tions . 

In hj s survey paper, S t -^tter [1978] exposes the gen-^ral 
structure of diff'"-‘rent tc^chniqu'^s for iterative 3.mprov ament 
based on I'l'^rr^'d correction or defect correction, and also 
prov'a th'-' convorgrnc 3 of iterates of the IlfeUmrthod towards a 
fixed point using contraction principle but avoiding a.symp- 
totic jxp-^nsions. 

As a continuation of Lindh'-rg's work SkGel[l980] 
disease's theoretical aspects of iterated deferred correc- 
tion avoiding th^ asymptotic ‘^xnansions of the viol'll 
error. 

In tae light of the past research outlined above, 
tho pressnt chrsis introduces ne-w and goncralio.cd methods 
of itoritDve improv Biriv^nt and presents a varient of IDeC as 
well as venous novel applications of IDeC . 

1,2 Suia .jary of the Thesis 

Th^ next chapter deals with now applications of IDeC 
based on quadrature methods for the solution of linear 
Fredholm, and linear and nonlinear Yolt-^rra int'^gral 
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equations of th^ second kind. e also bring out the 
difficulties encountered in th*^ application of IDpC to 
igcnv.alue problpms . 

Chapters 3 to 6 pi-osont four new acceleration t^ch- 
niquGS wliioli hav"‘ boon named as Successive Sxtra.polatGd 
II oC (oEIDeC), Iterated Deferred Defect Correction ( IDDeC ) , 
Successive Extrapolated IDDeC (SSIDDeC), end IDoC on Extra- 
polation (iDeCB). ThDse m^thodSy which, involve IDoC with 
thp US ; of interpolating polynomials, are described as 
apple "d to quadraturo methods of Fredholm integral equations 
Th ' 'ify^ptotic expansions for th'^ methods ar-^ derived* eenplo 
ui’^nt'^tion details are provided; ^nd illustrative numerical 
exa'iiples h*'ve been worked out. Applications to linear and 
nomine er Voltorra equations have also boon considprad. 

A v'^'^iant of IDoC call'>d Successive Updat'^d IDeC 
(SUIDfC) as appl L ;d uo nonlxnoar Volterra equations is 
pr'c^nt-^d j.r ChapuPr 7» A comparative study of all these 
acc '’I ' rat: on t ■'chniques is made in Gheptor 8. The ninth 
chppti^r IS d'jvored to applications of defect correction, 
with t-' " use of splines ^nd polynomials, to different types 
of pee’tiel differential equations. In the last chapter 
courluoc.ors are drawn and future research possibilities ore 
d 1 sous sed. 

The class of acceleration techniques has been en- 
riched remarkably by the IDeC method. When applied to 
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Vclterra equations , it generally produces solutions of high 
accuracy even for relatively large stepsiSies. For Fredholm 
equations with smooth kernels and functions^ the easy-to- 
if^plement SEiDeC method extends the order of convergence 
. bitrarily beyond the maximum order attainable by IDeCy 
'•ach iteration of iDeC and iDDeC increases the order of 
. jnvergence by n and n+p respectively ^ where n is the order 
ol the basic discretization scheme and p is the increase 
in order due to deferred correction. Because of its 
guaranteed convergence^ IDDeC is superior to IDC also. The 
r.iethod of SEIDDeC extends the order of convergence of IDDeC 

Numerical results confirm theoretical analysis and 
i.ndicate the capability of the defect correction techniques 
in providing general software for the automatic computation 
of solutions to problems not only in integral equations but 
also in other operator equations. 

The thesis has been written in the accepted mathematical 
style of establishing certain theoretical results and giving 
numerical confirmation by specific examples. It is well- 
known that most theoretical results are guided by plausible 
reasoning and numerical results indicating the plausibility. 
Chronologically/ the work presented in this thesis had gone 
through the plausible arguments, numerical computations and 
theoretical justifications. The details are given below. 



(i) Having found the success of Zadunaisky*s method 
in ordinary differential equations, it looked highly plau- 
sible that it should work in partial differential eq.uations 
as well. As such our initial work is what is contained in 
Chapter 9. However, due to serious limitations in our being 
able to create satisfactory bivariate function to compute 
the defect we could not complete this problem. So we 
attempted solutions by a series of univariate interpolatory 
functions but without much numerical success, 

(ii) ¥e then turned our attention to applications of 
IDcC to integral equations. We found that the use of 
interpolatory polynomials is preferable to that of inte3>- 
polatory cubic splines, because the former has a smaller 
interpolation error. 

(ili) The derivation of theoretical results as well as 
the construction of new acceleration techniq[ues described 
in this thesis were motivated by an analysis of numerical 
results obtained by IHeO methods for integral equations. 



CHAPTER 2 

APPLICATIONS OP IDeC TO IFTPGRAL EQUATIONS 


2. 1 Introduction 

In the area of application of the method of iterated 
dpf-^ct correction (iLeC) to ordinary dif f i=rnntial equations, 
substantial research has been done by Alfeld (1975), 
Zadunaieky (1976), the group led by Prank (1975-1977), 

1 tetter (1979), and Zadunnisky and Lafferriere (1980). 
However, th-^re have beeh only a few or no applications of 
iDeC to other types of operator equations such as partial 
differential equations and integral equations. 

In this chapter we deal with applications of the 
relatively n'^w ILeC method to int-^gral equations. We 
consider IDeC techniques Las'll on quadrature methods for the 
num'^rical solution of Fredholm and Volterra integral 
equations of the second kind. 

The next section contains a description of the 
method of ILeC as applied to linear Fredholm integral 
equations of tho second kind, and illustrates the effi- 
ciency of ILeC by means of numerical examples. Similarly 
Lections 2.3 and 2.4 deal respectively with non-linear and 
linear Volterra integral equations of the second kind. 

The difficulties involved in th^ application of ILeC to 
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eigenralU'^ prolleiUvS ar? brought out in Section 2.5. 

2.2 S redfaolm Integral E quations o f the Seco n d Kind 

Consid'^r a 3in'^ar Fredholrr int-^gral equation (FI^l) 
of th-^ second kind; 

a 

(2.2.1) y(x) -/k(x, t) y(t) dt = f(x), a <x < b. 

b 

In this equation th-' kernc'l k(x, t) is a known function, 
f(x) IS also known, and th^- function y(x) is to be detcr- 
nuned. y assume that k(x,t) is well-b<^hav^d for any x 
and t in [a,b] and f ( x) is sufficiently smooth so that a 
unique solution y(x) exists. Pn approximate" solution to 
y(x) can be found as describe>d below. 

let be a positive ev^n integer, h = (b-a)/]?, and 
x^ = t^ = a + ih, 1 = 0(1)F. he then use a Newton-Cotes 
quadrature formula such as Simpson’s rule to replace the 
integral in (2.2.1) by a sum. Th'= substitution x = x^ in 
the resulting equation leads to the equation 

N 

(2.2.2) - £ w^ ^ = 0(1)N, 

3=0 

where is an approximation to y(x^), and the w^ are the 
weights in the quadrature formula. The system of simul- 
taneous linear algsbraic equations (2.2.2) may then be 
solved for the determination of tb= Yj^ by G-aussian 
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elimination. The above method is called a quadrature method 
and is denoted by Quad(h). 

¥e improve these approximations by applying th^:* 
method of Iterated Defect Correction (iDeC) as explained in 
the next section. 

2 , 2.1 The Method of Iterated Defect Correction 

V/e first estimate the unknown global error 
y(x^) - Y^ in the numerical solution Y^. Assume that the 
integer N can be expressed as a product of two integers m 
and s, where m is a fixed even integer greater than two, 

¥e divide the interval [a,b] into subintervals 

II, I 2 , I 3 , where I^ = 3 -l)m’ ^jm^’ ^ = lil)s. 

Dor each 3 , l^t P^*^^(x) be the polynomial of degree m which 

t) 

interpolates the m + 1 points 1 

We then define for x in [a,b], the interpolating 
function as the collection of polynomials P ^*^2 

III J y ul 

such that 

P^°\ x) = P^^^( x) forxel., j = l(l)s. 

ITi J y 11] J 

(¥e observe that these polynomials depend also on the step- 
size h '^v-^n though the notation does not bf^ar out this 
fact. ) 
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¥-5 now construct a new problem 
b 

(2.2.3a) y(x) -j'k(x,t) y(t) dt = 

a 


whore 


a < X < b 


b 

(2.2.3b) g^°\x) =P^°^(x) -jk(x,t) P^°\t) dt. 

a 

SincG g^^^(x) - f(x) IS small, (2.2.3a) is a neigh- 
bouring problem (NP ) of the original problem (Op) given 
by (2,2.1). The exact solution of this ¥P is the known 

function P^*^\ x) . ¥e then solve this problem (2,2.3a) by 

the sam ^ method ¥uad(h), and donoto its solution at x^ by 

r] (0), 1 = 0(1)N. The error in this solution is known, and 
1 

IS equal to - h ^ • 


Since both Op (2.2.1) and NP (2.2.3a) are solved by 
the same method, we postulate that the unknown error 
y(x^) - IS approximately equal to th'’ known ">rror 

P^^*^\r.^) - Therefore, th-- approximation Y^^^ giv'^n 

by 

(2,2.4a) y[^^ = '^i + - "^ 1 *^^ 

IS n bett-^r approximation to y(x^) than Y^. This principle 
of jrror -estimation is due to Zadunaisky (1976). 
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The process above may be repeated as many times as 
are necessary. The improved approximations are then 

given by 

(2.2.4b) r = l,2,..., 


where the function x) , defined analogous to theP^*^\x) 

of (2.2.5a), interpolates the points (x^» Y^^^), i = 0(l)N 
( 1 * ) 

and the is the nximeriCal solution of the neighbouring 

problem 

^ (r/ 

(2,2.5a) y(x) -/k(x,t) y(t) dt = g(x), a< x< b 

a 


where 


b 

(2.2.5b) g^^^(x) =Pj^^^(x) -/k(x,t) P^^\t) dt 

a 

for which the exact solution is P^^^(x). This technique 

m 

IS called the method of Iterated Defect Correction 
(Sathiara^ and Sankar [1982] ) . 

( x* ) 

Ve observe that, for each r = 0,1,2,..., p^ is 
obtained as the solution of a new system of N+1 equations 

If 

(2.2.6) - £ w^k( x^, t^ )p^^^ = g^^\x^) 

0=0 


This system has the same matrix of coefficients as the 
system (2.2.2), but the term replaces the right 
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b.'^nd sid 7 i:(x^) of (2.2.2). Con s -q.u -ntly, wh.'^n finding the 
solution of (2.2.2), if r''prGsent the matrix of 
co-ef f ici-'nts in the product forr- LIT, where L and U ar® 
trnn'^ul'^r nitrico?, the Gaussian olinination n-^ed he 
lion^ only onc-^. This ro^uir^s 0(F’’) operations, whoroes 
th'’ dot '’rmim tion of th? solution of (2.2,6) req.uir3S 

only 0(11'’) operations, if the computation of the right-hand 
Side of (2.2.6) is ignored. 


Inpl ’'1 ntation 


hh 'n solving (2.2.6) for th'' numerical solution 
of tho n-'ighhouring problem (2.2.5a) with a fixed r, 
wo ne''d to compute the valu s g^''^\x^), i = 0(l)F. By the 

definition of P^^Hx), we have fror- (2.2.5b) 


(2.2.7) 


"b 

(r)/ \ J l(x.,t) P^^^(t) dt. 

g (x ) = Y" ' ^ i’ ^ m 

a 


The integrand of the integral in (2.2.7) is a known 
function, and therefore th'^ integral may conveniently be 
comput d as accurately ns n-=ccssary. In our computations, 
TT\G routines were used for this purpose. 


2.2.2 tum~-rical results 

In this section the relative performance of the 
m''thods of quadrature and iterated defect correction is 


illustrated. 
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Tho problems below were solved using 16-18 signifi- 
cant digits. The first four of these problems are taken 
from Netravali and de Pigueiredo (1974) and the last 
problem from Baker (1977). 

Problem 1 

1 

y(x) - j xt y(t) dt = - X 

0 

3xact solution; y(x) = e . 

Problein_ 2 

1 

y(x) - j xt y(t) dt = sin tcx - x/u 
0 

'^xact solution; y(x) = sin xx 
Probloa 3. 

1 

y( x) - j x'^ y(t) dt = X - x^ + x^ e^(l - x) 

0 

3xact solution; y(x) = x 
Problem 4 

1 

y( x) - J x'^ y( t) dt 
0 

exact solution: y(x) 

ProbloPi 5 

1 

y(x) - y(t) dt = 1 - ( - l)/x 

0 


xx'^(e^ + 1) 
= sin xx - 2 2 — 

( X + % ) 

= sin XX 


exact solution; y(x) = 1 
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Th? msthods were eppli^^d using both thf’ trapezoid 
and Siiipson’s rules. Ibo degree e of the interpolating 
polynonials was chos<^n to be 6, 8, 10 or 12, and the number 
F of intervals in [0,l] to be n, 2n and 4m. In Tables 2.2.1 
to 2.2.5 are listed the naximum values of th'^ absolute error 
in th 1 computed solution for th^ relevant cases. 

"/e denoti the ID'-iQ methods based on the trapezoid 
rule and Simpson's rule r ispmetively by IDeC^ and IDeOg. 

In our tables, thi results for the quadra tur'^ methods 
based on the trap Dzoid and Simpson's rule appear ns entries 
corre spending to the iteration number 0 of the corresponding 
IDoC n ,thods. An entry with asterisk ( *) indic'^tes that 
the subsequent it'^ratas also hiv the sa^m entry. 

The results for Probl<=ns 1 and 3 by IDeC^ with 
ui = 10 and 12 ar'i giv-n in Tabl'- 2.2.1, and those for 
probl'ins 1 to 4 by the IheOg in Table 2.2.2. Wo find that 
both the ID-iQ methods work much bett'ir than the underlying 
quadrature nithods and that IDnCg has a more rapid conv-ir- 
g'noi then IDi^Cg^. With a single defect correction, IDeCg 
IS capable of producing a highly accurate solution. 

In T'lbl'i 2.2.5, w(i display th-' r'^sults for the ID'iC^ 
and Ili'Cc. .’’'■'thods with m = 8, and in Table 2.2.4 those for 
the ni'ithods with m = 6. '-/a notice that thi iterates of 

both tl' . IDeC methods converge to the same solution, and 
that IDnCq needs generally just oii'i iteration to reach the 


point of convergence-. 



Table 2.2.1 Maximum errors -IDeC Method on Trapezoid Rule 
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Table 2.2.2 Maxi^mum Errors- IDeC Method Bas ed on Simpson' s Rale 
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Table 2,2.4 IDeC MethodsBased on Trapezoid (T) and Simoson's (3) Ruleswith 



8E-10 3,0E-10* 


Table 2«2«5 Maximum Errors in solution -bo Problem 5 by IDeC Method Based on Trapezoid Rule 








03 








tH 

tH 






1 

1 

1 


VO 





Pq 

CO 

9E 






• 

• 

• 






tH 


tH 




o 

03 


rH 

ro 



CO 

tH 

tH 

00 

tH 

tH 



1 

1 

1 

1 

I 

1 


LO 

w 

m 

w 

w 

pa 

pq 



03 

03 

tH 

03 

00 




« 

• 

• 

• 

• 

• 



VO 

tn 


tH 

cr\ 

O' 





o 



tH 



I> 

a\ 

tH 

t> 

c\ 

rH 



1 

1 

1 

1 

1 

1 



w 


w 

ra 

pq 

pq 



On 

ro 

m 

03 

o 

03 



« 

t 

• 

• 

• 

• 




00 

tH 

tH 

03 

ro 

u 








Q) 


VO 

l> 

Ch 

VO 

00 

(J\ 



1 

1 

1 

1 

1 

1 

!z: 

CO 

W 

w 

W 

(■a 

pq 

pq 


CTv 

on 

m 

03 

tH 

ro 



• 

• 

• 

• 

• 

• 

c 


CO 

tH 


tH 


rH 

0 








•H 








-P 








fO 








U 








0) 


LO 

VO 

r- 

m 


CO 

■P 


1 

1 

1 

1 

1 

1 

H 

03 

W 

W 


W 

pq 

pq 



rH 

03 


on 





• 

• 

IT 

« 

•r 

• 



m 

03 

tH 

! rH 

CO 

lO 




ID 

VO 


in 

VO 



1 1 

1 

1 

1 

1 

1 


rH 

i M 

w 

w 

pa 

pa 

pq 



m 

m 

in 

1 m 

r- 

CO 



1 • 


• 

• 

t 

• 



03 

m 


rH 

1 

tH 

03 



m 



i 

on 


un 



1 

1 

1 

1 

1 

1 


o 


W 

pa 

ca 

W 

S 




o 

in 

m 

00 



• 

• 

♦ 

• 

• 

* 



03 

VO 

tH 

tH 

fO 

CT\ 


2; 

00 

VO 

03 

O 

o 

O 



tH 

on 

tH 

03 



e 


CO 



O 








rH 



CO 

» 

CO 

CO 


o 

I 

w 

« 


I 

w 

o\ 

m 

if) 


in 

w 

m 



CN] 

tH 


O 

tH 

I 

w 


03 




CO 

I 

w 

VO 


in 


I 

w 

w 

CO 


VO 

I 







I 

w 

r- 

# 

03 


00 

03 ^ 


03 


6.7E-5 1.3E-6 2.5E-8 5.0E-10 9.9E-12 


23 


Tabl^ 2.2.5 giv'=>s tte results for Problem 5, an 
example for which IDeC works well ev'=n when there is a 
discontinuity in f(x) of ( 2 . 2 , 1 bei>^ c(^fn-ed d>7 1 

The t^st results also show that 1iie greater the 
degree m of the polynomials used in IDeC the higher the 
accuracy in its solution. 

Th'^orotical results on the asymptotic order of the 
IDeC^ and IDeCg methods for eq.uation (2.2.1) are given in 
Section 3*2.2 of th-^ ni^xt chapter. 

2 . 3 Monl in^^ar Volterra Integral 3quetions of the Second Kind 

Consider the following nonlinear Volterra int^^gral 
eciuetion (VI'C) of thc' second kind: 

X 

(2.3.1) y( x) - J k( X, t,y( t) ) dt = f ( x) , 0 < x <_ a, 

0 

where k(x, t,y*) is assumed to b-^ continuous in x and t, and 
to saticfy a uniform Iipschitz condition in y, and f is 
continuous. The ^q.uation has then a uniq.ue solution which 
IS d snot'^d by y(x) ( Davie [1962 ]) . 

An approximate nunorical solution can be found by 
th-"’ following procedure. Ve choose a positive e-v-^n integer 
K, and 1 )t h = a/l and x^ = t^ = rh, r = 0(1)11. let 
denote ail epproxioia tion to y(x^), and f^ stand for the 
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valui f(x^). In (2.3.1) we set x = rh, and roplac- the 
int'^gr‘^1 by e quadrature formula -virxth ax^propriately chosen 
weight •’unctior-. Th'n wo can determine from 

r 

(2.3.2) Y^ - h £ w k(rh,Dti,lj) = r = 1 , 2 ,..., 

3=0 ^ 

with Yy = 3 ^ 0 ) = fg. In general, (2.3.2) is a nonlinear 
equation in the single unknown Y^ and can be solved by an 
it^retive procedure. 

'''/ 3 have conoid ■^red three step-by-stop procedures 
tor tho nuii'''rical solution of th''^ Yolterra equation (2.3.1). 
They ar^ based on the trapezoid rul^, Simpson’s rule and 
nid-point rul'. 

Trapezoid rule 

In the simplest case, when the quadrature formula 
used IS the trapszoid rule, tho equations (2.3.2) have the 
f orru 

r 

(2.3.3) ^r ^ k(rh,3h,Y^) - ik(rh,0,YQ) 

3=0 

- ^k( rh,rh,Y^), r =r 1,2,... . 

Starting with Yq = fg, wo solve the first of the 
equations (2.3.3) for Y^, and then solve th^ second 
equation for Y 2 a-nd so on. 
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IToble's method 

Por more accurate computation we use higher order 
Quadrature formulae such as Simpson's rule. In Uohle's 
method (Noble [1969 ]) , we begin with two starting values Yq 
and Y-]^ and compute, for 

r = l(l)N/2, from 

(2.5.4a) 


h r 
f _ + - Z [k 


Oy> O*? 1 


2r 2r ^ ^^^L^2r,2i-2 ^^2r, 2i-l ^2r,2i 


and Yr. , from 
2r+l 


.3-^b) 


h 


r-1 


^2r+l " ^2r+l ^ ^ t^^2r+l,2i-2 ^^2r+l,2i-l ^2r+l,2i^ 

^ 1=1 

5h 

^ ““ Lk2y^l^ 2r-2 ^ ^^'2r+l,2r-l ^ ^^2r+l,2r ^ ^2r+l,2r+l^ 


with 2r+l<N and k - = k( x , t ,Y ). 

^ D 0 

Eepoated Simpson's rule is used in (2.5.4a), whereas 
repeated Simpson's rule supplemented by the 5/8ths rule in 
(2.5.4b). The stability of the above method was proved by 
Noble [1969]. 
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Mid-point rule 

¥e may also compute the approximate values for 
y(Xj^) ^y the use of the mid-point rule as below, 

^0 = ^0 

■^1 = fi + hk(x^,XQ,YQ) 

Y 2 = ^2 ^ 2hlc( X 2 , x^, Y^) 


(2*5.5) 


"^2 1-1 ^ ^2 1-1 ^2i-l» 

1-1 

+ 2h Z 
3=1 

1 




i ^2 i'^23-1’Y23_i) 

3=1 


This is an explicit method. 

2 . 3.1 Application of IDeC 

The application of iterated defect correction to 
Volterra integral eq. nations is similar to its application 
to Fredholm eiuations in Section 2.2. The details are 
essentially the same, 

Ve first solve the original problem (2,3.1) by a 
specific method M(h), and find the discrete solution 
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Y^, 1 = 0(l)lT. As in Section 2,2, a smooth global approxi- 
mation \ x) is defined In [0,a], by extending the 
discrete solution at x^^^ in [0,a], We then construct the 
following neighbouring problem for (2,5,1): 

X 

(2.3.2a) y(x) - /k(x,t,y(t)) dt = g^^\x), 0 < x < a, 

0 

where 

X 

(2.3.2b) g^°\x) =P^^^(x) dt. 

0 


This problf^m is solved by the same method M(h), and the 
sclut-uOn produced is denoted by i = 0(l)W. We then 

compute a better approximation to the true solution of 
(2„3*-''-) by using the equations 


= Y. 




0(1)W. 


The nrocess above may be repeated, and the solutions 

( 2 ) ( ') 

Y^ , 7^. ... may be found as in Section 2.2. 

betails about the implementation of the method of 
Il)f^C on quadrature methods for the numf=rical solution of 

(2.3.1) are similar to those in the case of Predholm equation 

( 2 . 2 . 1 ) . 



28 


2.3.2 Numerical results 

The results of some computational experiments are 
given in Tables 2 . 3.1 to 2.3.9. "the computations were 

done using 16-18 significant digits on the DBG-IO 9 O computer 
at I.I.T. Kanpur. We programmed the three quadrature 
methods based on the formulae (2.3.3), (2.3.4) and (2.3.5), 
and the IDeC m«^thods on these three quadrature methods. 

We Used each method to solve the following five problems 
defined for x in [0,l]. Problems 1 and 3 are taken from 
Netravalifl973], Problem 2 and 5 from Grarey[l975] , and 
Problem 4 is from Hock [1981]. 

P roblem 1 

X 

y( x) = f X - x^ ) - x^ + x^ +j x^ y(t) dt 

0 

exact solution; y(x) = x 
Problem 2 

/ N X , x/2 

yix) = 3 - x+xe 

exact solution: y(x) 

Proble m ^ 

X 

— X 2 

y(x) = e” ( 1 + X + X ) - X +j’ xt y(t) dt 

0 

-X 


X 


xt y^( t) dt 


0 


= e 


X' 


exant solution: y(x) z= o' 
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ProTslem 4 

X 

y(x) = e“^ +j e'^"^(y(t) 4- dt 

0 

exact solution: y(x) = ln( x + e) 

Pi *o'b~' e m 5 

y(x) = 2 - 5x + (x + l)(5x - 1) 

X 

+ 2( X + l)(l - x) log(x + l) - 4 j'log(x - t + 1) y(t) dt 

0 

exact solution; y(x) =1 - x. 

Proliem 1 to 5 were solved for N = m, 2iii and 4m with 
Cl ^ 4« 6, 8,, 10 and 12. In Tables 2.5.1 to 2.5.9, we have 
listed rid'^: each method the maximum values of the absolute 
erro'^ between the exact solution and the computed solution 
for different problems for th'^ relevant values of m and N. 

IDoC methods based on the midpoint rule (m), trapezoid rule 
(T), and Foble's method (N) are denoted by IDeCj^, IPeO^ 
and IDoCj^ resnectively . 

I'h'^ errors listed in Tables 2.5.1 to 2.5.5 show 
that IDeC produces very good results in each of the three 
diff''-'‘^nt quadrature schemes. 

The relative performances of IPeCjyj, IDeC^ and IDeC^j. 
are displayed in Tables 2,5.4 and 2.5.5. The rate of 
increase in accuracy per iteration is rather slow for 
IDeCjj when compared with the rat^s for the other IDeCs. 



T^lo 2.3.1 IDeC Method Based on Trapezoid Rule 
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Table 2.3.6 IDeC method with m = 4 and h = 1/m 



Table 2.3,8 Modified increment methods, and the IDeC^ 
method with m = 10 


Problem 

N 

Method 

A B 

IDeC^ interation number 
12 4 

6 

2 

10 

1.9E-7 

2.0E-7 

1.6E-6 

3.3E-9 

1.9E-9 

1.9E-9 


20 

1 . 6E-8 

1.7E-8 

l.OE-7 

2.3E-11 

1.8E-12 

1.9E-12 

5 

10 

1.6E-7 

6.2E-8 

2.9E-5 

5.1E-7 

1.4E-10 

3.7E-14 


20 

8.4E-9 

7.0E-9 

1.9E-6 

8.2E-9 

1.7E--12 

6. IE- 14 


Table 2.3.9 

Spline method of Netravali, and the IDeC^ 

method 


with m = 10 


Problem N 

, IDeC„ Iteration number 

Netravali T 

12 4 

6 


10 

8.8E-5 

1.5E-4 

7.1E-6 

2.6E-8 

1.2E-10 

20 

8.8E-7 

9.1E-6 

l.lE-7 

3.4E-11 

1.7E-12 

10 

l.OE-5 

3.8E-6 

2.4E-8 

4.6E-12 

3.5E-12 

20 

8.7E-8 

2.3E-7 

3 .OE-lO 

1.3E-12 

1.3t-12 


2 
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is not a marked difference Toet-ween the performances 
of IDeC^ and IDeC^^^ except that to obtain a given accuracy 
the number of iterations requir'd by is smaller than 

the number required bv IDeCQ,. 

Table 2.3.6 illustrates the efficiency of the IDeCj^ 
method for a relatively large step size, and Table 2.3.7 
that of the IDeCgi and IDeCp|. methods -with the use of a 
relatively larg^ degree of polynomials. 

The results obtained by the IDeC^ method based on 
the trapezoid rule are contrasted with those by the modi- 
fied increiflont methods A and B of G-arey[1975] in. Table 2.3.8, 
and with those by th'^ splin'^ approximation method of Fetra- 
vali[l973] m Table 2.3.9. The ICeCg, method produces much 
more accurate results than any of the basic discretization 
methods even for relatively small values of F. 

■'/e conclude this section with th'’ observation that 
th'’ performance of IDeC methods based on quadrature schemes 
for non-linear VIEs is indeed very good in practice. 

2 . 4 li near Volterre Integral Equations of the Second Kind 

Fe now consider the numerical solution of the 
following Volterra integral equation of the second kind: 

X 

(2.4.1) y(x) - /k(x,t') y(t) dt = f(x), 0 < x < a. 

0 
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¥e assume that kCx, t) is continuous for 0 < x,t <_ a and 
that f(x) is continuous for 0<^ x < a. Then a continuous 
solution y( x) exists (Baker [1977] ). 

'll 3 choose a positive even integer F and let h = a/N. 
For r = 0(1)F, we substitute x = rh in (2.4.1) and use a 
q.uadrature formula to approximate the integral 

rh 

j’k(rh,t) y(t) dt. 

0 

This leads to the equations 

r 

(2*4.2) Y(rh) - 1 w , k(rh, 3 h) T( oh) = f(rh), r = l(l)N, 

0=0 

where T<' r*h) denotes the a-op'^cxime t: on to y(rh) and th'’ w 

r y J 

are th,' weights escociei^'^u wj bh quadrature formula used. 

Thus, starting wj I’l i.’iw initial approximation 
T(0) - y(0) = f(0), vf''' obtain ch'^ equaxions 

- w^ ^ k(b,0) Y(0) + k(h,h)] Y(h) = f(h), 

aoooe<**i>*ott« 

( 2. 4.3) 

1' -1 

- g W r , k(Hh, j'.o Y(]h) 

>'=0 

+ [1 _ k(Fh,Fh)] Y(Fh) - f(Nh) 



I'rom tile above successive eq.uationSj, the values 
Y(h), Y(2h), Y(Nh) arc calculable one after another. 

We will call this a uuadrature irn-'thod^ since we have 
adapt-^d here th"^ idea of th'^ quadrature method for Yred— 
holm oquauions. Mor'^ details on (2.4.1) are given in 
Baher[l977]. 


A-p-plico t3 on of IDeO and Fu rical results 


Th^ details regarding the application of IDpG on 
the quadrature m'^thod hosed on the trapezoid rule for th'^ 
linear VIE (2.4.1) in Section 2.3.1. except that now the 
sth neighbouring probdoirt "or (2.1.1) is of the form; 

X 


(2.4. 4) 


y(x) -,'*]:(x, b) y(i-) dt 





yi. 



dt. 


0 < X < a. 


W- .mplernented the o^uadratui"' and IDeC methods 
applicable to (2.4.1). We giio here results from the 
computational tests on the following sample problem taken 
from Baker [1977]; 

Problem 

X 

y(x) - 2/ cos(x - t) y(t) dt = 
c 

xnct solution; y(x) = 9^(1 + x)^. 
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More results on a few other sample problems are included in 
Chapter 8. The calculations were performed on the ICL 190 4S 
computer at Banaras Hindu University, Varanasi in single 
precision (9-11 decimals) arithmetic. 

¥e choose the number H of intervals in [0,l] to be 
4(4)16, and the degree m of the polynomial elements of the 

function ^( x) in (2,4.4) to bo H. Maximum absolute 

ei*rors in the computed solution by th® quadrature and IDeO 
methods are given in Table 2.4.1 below. 

Th^ first it-^rat'^ of the IDeC method based on the 
trapezoid rule is as good as the result obtained by the 
m'^thod which us.^s the weights of repeated Simpson's rule 
with th^ 3/8ths rul^^. The latt-'^r method gives the errors 
3 X 10 ^ and 2 xlO ^forU = 8 and S’ =16 respectively 
(see pag3 785 of Baki^r [1977] ) . 

Table 2.4.1 The J.oad rat ure a nd IDeC Qv-^thods bassd on tra-po - 
zoid rule 


‘-Quadrature 


IS'^0 method iterations 


1' 

method 

i H 

I 

i 

2 

3 

4 

5 

4 

4 . 4E-1 

5.8B-2 

3. 6E-3 

9.9B-4 

7. 6B-4 

7.93-5 

8 

l.lU-1 

5.0E-3 

2. OB-4 

2.1B-5 

2. 43-6 

2.73-7 

12 

5.1E-2 

7.1B-4 

4.6B~5 

4. 8E-6 

5.13-7 

5. 53-8 

16 

2. 9E-2 

2.5B-4 

5. 53-6 

4.13-6 

1. 23-6 

4.33-7 



41 


1 ) "'rf orinanco of IDeC on linoar Voltorra integral 
equation (2,4.1'' is quit.'- "ood, even for relatively large 
values of h. 

5 ’liic'nvalu’^ Pr o_'bl ^riS 

Ln this s'^ction w ? consid'^r the nunerical solution 
of ti:.' FrecHioln integral equation of the third kind given 
by 

b 

^2„5,1) ' k( x, t ) y( t ) dt = ly(s), 4 i b, 

where tht k rnel k(r,t) is given, and the aigenv-^lues 4 
and th^ corr-^spording '^ig'^nf uncxions y(x) ari^^ to be cal- 
cula bod, 

AS in Section 2,2, a sniple nunerical m'^thod for 

(2. 5.1) IS to repl'ce bh'^ inb'''gral by a quadrature rule, 
and collocn.te it th*' points x = x ^ . This leads to 

t: 

(2.5.2) 2 w_ k(x^,t )Y = A T i = 0(1)F, 

3=0 

where Y and A stand for an approximation to y( x^. ) and 

X X 

r ■'Spec tively. In matrix nototion, this b'-com'=s the matrix 
-ig'nv.ilu- problem 

(2.5.3) KY = AY 

Ip 

K = ^^"^3 ] tXiid. \ — o # » ^ ^ • 
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Solving (2.5.5) we obtain U+l eigenvalues , r = 0( l)rr 
and the corresponding eig^^nvectors Y^. The above method 
IS the ‘quadrature method’ for the eigenvalue problem 

(2.5.1) (Baker[1977]). 

A method for iterative improvement 

Since the application of the method of IDeC to 
eigenvalue problems is not straightforward, we shall now 
describe a modified method for estimating the errors and 
improving the results for the eigenvalue problem (2.5.1). 

Suppose the quadrature method based on the quadratin' e 
rule fi- produces an •eigenvalue of the integral, equation 

(2.5.1) , and th-> corrospording eig'^nv^ctor Y^^ ^ which 
approximates th.^ r^quir<^d '"^-igenf unction, let the function 
P^°\x) b'^ a smooth global approximation to the eigen- 
function corr'^sponding to th'’ eig'^nvalue ^ of (2.5.1). 

This function is computea by fitting a polynomial to the 

data points ^ least square sense. Then 

consider the problem 
b 

(2.5.4a) j’k(x,t) y(t) dt = -A^yCx) + g^^^^x), 
a 

wh^re 

b 

(2.5.4b) g^'^^x) =j k(x,t) t) dt - i^P^°\x). 

a 
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The exact solution for the above problem is P^^^x). 
notice that (P,5.4a) is not an eigenvalue problem of the 
form (2.5.1), but is a linear Fredholm equation of the 
form ( 2. 2.1) . 

¥e solve (2.5,4) by the corresponding quadrat\ire 
method based on the quadrature rule R, and denote its 
solution by Adapting the idea of Zadunaisky, ve 

expect given by 

(2.5.5) 


to be a better approximation to the eigenfunction than 

. The use of (2.5.3) will then give an improved eigenvalue 

,( 1 ) 


This process may b-> repeated for the determination f 


of (T 


( 2 ) ,\( 2 ) /^( 3 ) ^ 13 ) 






One difficulty in this method is that the integral 
equation (2.5.4a) is ill-conditioned, because of the 
eigenvalue of K. (See Noble [1977], p 919). 


2. 6 C onclusion 

We h'^ve successfully applied the m^^thod of iterat'^d 
defect correction on the quadrature methods for linear 
Fredholm, and linear and nonlinear Volterra integral equa- 
tions of th= second kind. 
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In the case of Predholm equations, our computational 
experiments show that the c^^eator the order of the basic 
quadrature method, the faster the converg^^nce of the IDeC 
method, and that an increase in the degree m of the poly- 
nomials yields an increase in the accuracy of the IDeC 
solution. The properties exhibited Iby our nimierical results 
paVe the way for our theoretical analysis of the iDeC method 
presented in the next chapter. 

As for the application to nonlinear Volterra integral 
equations, we have considered IDeC methods based on quadra- 
ture methods which use the mid-point rule and the trapezoid 
rule and the IDeC method based on Noble's method. IT wo 
arrang ' th'^so throe IDeC methods in th'^ increasing order 
of accuracy in their solutions, they will appear in the 
sequenc''- m-ntionod abov<-„ However, the IDeC method which 
uses the trepezoid rule is th^ niost preferable of th^ three 
methods since its underlying quadrature e '^thod is self- 
starting, It is remarkable that the IDeC methods fbt'n " 
given step-size generally produce more accure-te solutions 
than any other discretization m''thods with th"' same step- 
size. 

Again in the case of linear YIDs, the IDeC method 
based on th^ quadrature method that uses the trapezoid rule 
yields highly accurate solutions. 
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However, application of the method of IDeG to 
eigenvalue problems is not straighiforvfsrd and presents 
difficulties in contrast to Richardson extrapolation and 
deferred correction techniques which are easily applicable. 

On the whole, IDeG methods are easy to implem'^nt and 
very effective in practice. Moreover, they produce highly 
accurate solutions to problems in Fredholm and Volterra 
integral equations of the second Icind. 



CHAP TER 5 


THE H'jlTHQP OF SIJGCESSIYE EIITRAPOLATED IT3RATSD 
PEREC T CORRSCTIOH 


3 i In troduction 

In Section 2.2 of the last chapter ue described an 
application of iterated defect correction on the (Quadrature 
methods for the numerical solution of Rredholm's integral 
equations of th-^ second kind. In this chapter we derive an 
assmptotic expansion in even powers of stepsize h for the 
global error in the solution produced by the IDeC nathod, 
and show how the numerical results of Section 2.2.2 confirm 
our theoretical analysis. 

Applying Richardson extrapolation repeatedly on the 
IPeC method, we present in Section 3.3 the technique of 
Successive Extrapolated Iterated Pofect Corri^ction (SEIP^C) 
r. ad the resulting asymptotic expansion for the global 
->rror, Eum-'^rical examples are given in Section 3.4 to 
d^ronstrate the superiority of SEIDeC over IDeC. 

The implementation of SEIDoC is easy and gives a very 
favSt convergence for smooth kern. Is and functions. 
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5*2 As^yiripto tic ~^'Xpansion f o r ID eC 

An application of tL^ principl'= on the quadrature 

methods for th' second kind ^ linear Fredholm integral equa- 
bicns (2„2„1) was discussed in S'ction 2.2 of the previous 
cliap'c'T'. Iiow an asymptotic expansion for the global error 
in the soli Lion produced by the IDeC m‘'^thod is derived in 
this section. 

l^t g(x^) = and donot'^ the column vector 

m _ 

( P’.Qf ” ’ ' ’ th-^n ropr'^sent the linear 

svstem of equations (2.2.2) 

iv 

C=0 

in matrix notation as 

(7.2.1) (I - K)T = f 

wh'^re I and K ar^ (F + l) y (P + 1) matrices. But the voc- 

T 

tor Y = ( y( xq ) , , . . , y( Xj^t ) ) satisfies 

f ;.2.2) (I - K)Y = ± + r 

^jher^ th'-- -^rror x in the quadrature'' formula is given by 
b _lv 

(7.2.7) T = lk(x.,t) y(t) dt w k( x t ) y( t ) . 

d. J ±. J J. J J 

a 3=0 

Substracting (7.2.1) from (7.2.2), we have 


(I - K)(y - T) = T 
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^!ow il tht^ non- singularity of (l - Pi) is assumed, th^ ^^rror 
in th.e computod solution can "b ' -writton as 


(3.2.4) 


y - T = (I - K)"^ T. 


similarly, for any fix'll r = 0, 1, 2, tiio error in 

uho computod solution 


ri 


(r) 





T 

) 


for Np 

(3.2.5) 

( 3 . 2 . 6 ) 


(2.2.5) or (2.2.5) is 


'5(r) 


' m 


- p 


(r) 


.(r) 


Ij 


= (k(x „t) P 


( r 


J 

a 


m 


given by 

(l ~ Y) ^ whero 

F 

it) at - 

3=0 


Fow for the quadrature method Quad^(h) based on the 
trapezoid rul-', th. weights w in (2.2.2) are chosen such 

J 

that 


Wq = Wj;^ = h/ 2 . w = h for j = l(l)F-l, 

u 

and the quadrature error in (3.2.5) takes tho form (see 
Play 'O'rs [197^ ] walsPi [1974 ] ) 

(3.2.7) = A2(x^)h^ + A^(x^)h^ + Ag(x^)h^ + 


7o assume xhat (l - K) is stable in th^ ssnso of 


Foble (see Nobl'^[l977] , pp. 921-922). Th“n bpcauss of 
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equations (3,2„4) and (3.2.7), th'^ order of the error 
7^ - y( x^) in the approxnrate solution is the same as 
the ordjr of the quadi^aturo error, and we can write 

(3.2.8) - y(x^) = B2(x^)h^ + B^(x^)h'^ + Bg(x^)h^ + ..., 


L't Both Y(x) and Y(x,h) represent a continuous 
approy ination to y(x) of (2.2.1) constructed from the 
discr'^'te solution Y^ obtained by the method '4uad^(h). 

Th'^n, because of (3.2.8), wo hare 

(3-2.9) Y(x.h) - y(x) = B 2 (x)h^ + B^(x)h'^ + Bg(x)h^ + ..., 
whf’re B2 (x), B^(x), etc. are suitable functions. 


Now, for each 3, the expression P^^\x) - Y(x) is the 

iHy 3 

error in the polynomial P^^\x) duo to int'^rpolation (see 
Ralston and Rabinowitz [1978] ) of the points ( x^, Y^) with 
X e I (the moaning of I as in S'^ction 2.2.1). Since 
P^^(x) = (B^j(x) - Y(x)) + Y(x), wo have from (3.2.9), 


( 3 . 2 . 10 ) = y(x) + B 2 (x)h^ + 0(h'^). 


Consequently, for the BP (2.2.3), the equation (3.2.6) 

yields 

b £_ 

= T + h^ [ jk(x., t)B^( t)dt - w k( X , t )B2( t ) ] + O(h^) 

3=0 

= T + O(h^), with tha aid of (3.2.7). 
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.( 0 ) 


\! 3 can therefore write ' as 


(3.2,11) = A.(x )h^ + a.(x )h^ + a.(x )h^ + 

1 4 1 D 1 


Moreover, from (2.2.1), (2,2.3To) and (3.2,10) we have 


^0), 


"^^^1^ “ ^^^1^ = 0(h2). 


It follows from (3.2.5) that 

( 3 . 2 . 12 ) ^ 

which IS analogous to th^ expansion (3.2.8) for - y(x^), 


The us-^ of Zadunaisky's principle 


r(l) 


^r' = - I:. 


( 0 ) 


yields th'^ equation 

(3.2.13) y[^^ - y(x^) = B^^)(x^)h^ + B^g^\x^)h^ + ..., 

with = B - p 3 = 4, 6, ... 

In a similar manner we can show that the following 
resu3+ holds, assuming that both Y^ nx,h) and Y(x) denote 
th? approximation to y(x) of (2.2.1) produced aft^r the 
r oh iteration of IDeC. 

Bor r = 1, 2, ..., we have 

(3.2.14) Y^^\x,h) - y(x) = B^^\x)h'^ + B^^^vx)h'^^^ + 


0 • «i ^ 
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and 

^3.2„15) - f(r^) = Oda'^), 

where q = aiin( 2r +2, ra + 2) ; sinc^ for any fixed r, we 

hav*-' 

3 = 0(1)U, 

and fn.rofor. (see Ralston and Rabinowity [1978 J ) 

b 

j k(x^,t)(P^" \ t) - Y^^\t))dt = 0(h“^^^). 


Thus in the case of the trapesoid rule, oach iteration 
r, with 1 T < m/2, of th'^ IfeG increases th'^ ord'^r of 
accuracy in the solution by two, and th= asymptotic expan- 
sion for the global error in the final solution 
of the IfcG method is given by 


( 3 . 2 . 16 ) 


Y ( 
ITeC 


y( r) 


Y+7 


Tf the quadrature method 4uad( h) is based on Simpson's 
rul3 instead of the trapezoid rule, the t -rm containing h^ 
does not appear in eo^ nations (3.2.7) to (3.2.10), and the 
relation b ^tween the errors ir the approximate solution of 
th ' original problem and th^ first neighbouring problem 
encountered in the defect correction procedure is given by 


( 3 . 2 . 17 ) 



+ 0(h® ) . 




u: 
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Consoq.u ■'nily, wo can show that 

(3.2.13) - y(z^) = + ... 

and that the resulting asymptotic expansion for th' ID'=0 
method remains th-^ satii= as (3.2.18). 

In fact, for a fixed la, th"' asymptotic expansion 
(3.2.16) for tho global error in th-' solution of (2.2.1) 
is indop -indent of th'-^ underlying Howton-Cotf^s formula which 
the ID''G m-'thod uses. Thus tho "rror is of ord-'r h'^"*'^, 
where h is the st-'psize and m is th."^ dogreo of tho poly- 
nomial employed in th'^ n'-'ighbouring problems. 

Thus v(? have th'^ following th-^orem: 


The orem 2 . 1 

If tho order n of th' quadrature method for (2.2.1) 
IS 2 or 4, and the dogrcB m of interpolating polynomials 

( I* ) 

IS oven, then the solution of the rth iterate of ID-^C 

satisf I'^s 


( 3.2.19: 


r( r ) 
X 


y(x^) = 0(h 


min(n + rn, m + 2) 


r = 0, 1, 


/n+l -1 


Iloreovor, th'-^ asymptotic expansion for the global orror in 
tho final solution of Yjp^q(x, h) of th'-^ ID^-’C has the form 
( 3 . 2 . 16 ). 
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uia.^r ic a l r>?sults 

For th ^ rosults on ord''r of conv^rgGnc-^ for th'^ 
ILoC n'tliod applied to th'^ t '^st-protlens in Frodholm 9 q.ua- 
tionOj w- rof^r to tho tablos in Section 2.2.2, and comput'^ 
tlic rauios of th > actual computed errors in the solution 
by th Ji)?0(i,h) ano IDoCii, -^h) m^^thods. Th''se ratios 
d -"noted Dy V(i,h) arc d^fir-^d as follows; 

•^rror by ID^C(i,h) 

(3,2.20) ^(i,h) == T-- , 

orror by ID8C(i,-^) 

\7h^"r--' a xs th"- it -'ration nurib^r ard h th-" steps ize. From 
th ■'SO vilu 's th-' orders of IDeC are calculated for a f^w 
oxampl'-s. Table 3,2.1 lists these values and the ord-^rs 
of ID'iCs based on th' trapozoid j’ulo, and Table 3.2.2 
those las')! on Simpson's rule. In both th" tablos L 
denol-^s th " face that th^ machine accuracy has bo-^n reached 
for at li ast on- of th" Ix/oC solution value's from which 
I'Cioh) has b-"'"n calcu3.ated 

Th- T-'-sults ir Tables 3-2.1 and 3-2.2 on th" orders 
of conv"rg"nce for th'" IDeC iterates ar'" in accordance with 
Ob'" thcoriGical results of Thoor'm 3.2.1. 



Tab 


3.2.1 The values V(i,h) of (3.2.20) and the orders 


ror the ich iteration of ijOeC^. 







Iteration Number 



i'L' 

Problem 

h 

C 

1 

2 

3 

4 

5 

6 



1/3 

n 

16 

65 

253 

997 

997 


8 

1 

1/6 

4 

16 

64 

255 

1009 

942 




C rder 

1 

4 

6 

8 

10 

10 




i/iO 

4 

17 

67 

273 

916 

3296 

3207 

10 

4 

1/20 

U 

16 

65 

260 

978 

4120 

4134 



Order 

2 

4 

o 

8 

10 

12 

12 



1/1? 

4 

16 

64 

231 

16566 



12 

2 

1/24 

4 

16 

64 

257 

L 




Order 2 4 6 


8 


14 
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b l - 3.2,2 Th" valgg^.s Y(i, h ) and th^ orders for 
th^ Lth -teration of IDeOg„ 


m P3’0bl-’m 

h 

0 

It. rati 

1 

.on Fumbi^r 

2 3 

o 

l/lO 

16 

236 

3217 

5213 

1/20 

16 

254 

4156 

4164 


Order 

4 

8 

12 

12 


1/12 

16 

250 

40 23 

L 

C\J 

i — 1 

1/24 

16 

254 

L 

L 


Ord or 

4 

8 

12 



Howovcr, ^n th. Craso of Problom 3r comput'^d 

oxdors for IDoCij with, rn = 8 wor. found to h'-' 2,4,6,8,10, 
r'liid 12, ond in cas-^ of Problom 5 thoy wer-^ 2,3,4, 5,6,7 and 
0 ( rr-f r Soctaon 2,2.2). Tho extro ord ^r of conT-^rgonce 
in tho form'''! caso is attributed to the fact that since the 
pxact solution is y(x) = x, the polynomials P^^ which 
ort^nd th-- discrot"' solution do not introduce an error due 
to int ■'rnola.tion. In thi^ latter case, thi^ increas'^ in th-^ 
order of conwergenco is not two, b'^causo of the disconti- 
nuity in tho function f(x). 
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5* 5 Suc cessive Extrapola t ed Iterated Defect Correc tion 
(SaiDeO) 


3,3.1 Asy mptotic ex-par-si ori fo r B^ iPeC 

low that th°re exists an asymptotic sspansion (3.2.16) 
tor the ID“0 metbod on the quadrature methods for the linear 
ii'redholin equation (2.2,lj, Richardvson extrapolation or 
'deterr-^'d approach to the limit' (see Ralston and Rahino- 
-^_tz[l978]) may be used as explained below. 

By computing the approxiraaxion 'by bhe IBeC 

method for different values of h^ and combining two approxi- 
mations which have similar errors, we can eliminate succe- 
ssive error berms 3=2, A, ... in (3.2.16), and thus 

obtain more and more accurate approximations. 


Bor instance, suppose we obtain an approximation 
^11) oO (r,h) with a particular stepsize h, and then halve the 
stepsize and get a second approximation yb), and 

1 inally find a new approximation 7'jpgQ(x,h) by combining 
the first two approximations as in the equation 


.,m+ 2 


(5-5.1) 




f-n', _ri(x, sb) — Yjp^Q(x,h) 


•IDeC' 


~,m+2 


Then w^ obtain 

(5.5.2) Y^j^gQ(x,h) = y(x) + G^_^^(x)h®^^ + 0(h“^^), 
which has no h^"^^ term. 
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If we now determine a third approximation, again by 
halving the stepsize -^h, w^ can combine it with the second 
to obtain another new approximation x, -jh) containing 

an error expressible as + 0(h'^'*'^). Prom th = 

new valui^-s x, ^h) , we can combine a 

still anoth-r new approximation Tjx)^^(x,h), from an equation 

similar to ( 3 . 3 . 1 ) but with th-^ exponent m +4 instead of m+ 2 . 

2 

ITow the -expansion for the approximation Yjp^Q(x,h) has no 
C^+ 4 (x)h'^^^ term. 


In general, after the elimination of the first g 
coefficients in (3.2.16), the form of the as 3 nmptotic 
expansion is 

(3.3.3) ^ 




Thus the principl'^ of Richardson extrapolation may 
be pcrfonn'^d successively, using the solutions produced 
by the IDeC laethod for various stepsizes h. Ive call this 
technique the method of Successive Extrapolaxed Iterated 
Defect Correction (SEIDeC). 


3 . 3.2 Implv’ ^m entation 

The method SSIDeC is convenient and efficient for 
automatic computation of the numerical solution of (2.2.1), 
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provided the unknown solution y(t) and th“ k^rnol k(x,t) 
ar-' sufficiently smooth. 

Calling- the process of computing tho values 
as th? .ith order 3xtra-Dolatio n, we illustrate the procedure 
foj computing tho SSIDeC solution for (2.2.1) up to an extra- 
polation of a given order until the valu'^ from the 

lat st extrapolation compl-^ted differs from the most 
accurat-^ value of previous extrapolation hy less than a 
tolerance TOL. 

Begin by calculating the IDeC solution values for an 

h 

appropriate stepsize h and then for and entering them 
along with the solution from the first order extrapolation 
as in thf^ following table; 


S tepsizG 
Solution 
1st order 


4-h 


■IDeC 


(x,h) 


•IDeC 




'xtrapolation 




If the condi'.on j ^ 

satisfied or if = 1, then take Y^^^^(x,h) as the 

r-'q.uirod value and t<=rminate the procedure. Otherwise, 

ll 

find the IDeC solution with stepsize and expand the 
table to th-^ following; 
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'• 1 ' ' 1 

Stepsize h. jh. -^h 

Solution "^IDeO^ ^ ^IDeO^ "^IDeC^ 

1st ord'^'r extrapolation ^Tnori^ 


2nd order extrapolation 


"IDeO^ 


If either of the two conditions 

~ ^ ^max “ ^ is satisfied, 

2 

take Yjj^^^Cxjh) as the output and terminate the procedure. 

Otherwise find the IDeG solution with stepsize -g, and 
expand the above table again. 

In the computer memory this table is represented as 
a two-dimensional array of the form 


> 

o 

o 

o 

H 

> 

O 

ro 

■^10 

^11 ^12 

^20 

^21 

^30 


where A. 

J- j 

^IDeG^ 


3. 4 Hum eric al Results 

In this section we give computational results on the 
application of SEIDeG on q.uadrature methods for Fredholm 
integral eg.uations of the second kind. We consider here 
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Problems 1 and 2 of Section 2,2.2. 

Th'' degroo m of the polynomials used was taken to be 
4 and 6. Maziimum absolute errors in the final solution 
produced by the IDeC and SPlIDeC methods with m = 4 and 
m = 6 ar given in Tables 3*4.1 and 3.4.2 respectively. 

In these tables the SEIPeC method with the ith order 
extrapolation is denoted by SEIDeC(i). 


Tabl-' 3.4.1 IleO and SEIPeO methods with m = 4 


e,i 


12 


Method 


Problem 1 
24 


48 


Problem 2 
12 24 


48 


IDeC 

SEIDcCd) 


SEIDeO( 2)9. 5S-18 


il,2B-8 1.9S-10 2.93-12 3.43-7 5.13-9 8.03-11 

6.83-13 2.73-15 , ' 1.43-10 5.53-13 


1. 63-14 


Tabl 

- 3.4.2 IDeO and 

SBIDeC 

methods with m = 6 


N 

M ethod 

Problem 1 

12 24 

48 

Problem 2 

12 24 48 

( 

iDeO 

5.53-11 2.13-13 

8.23-16 

1.23-8 4.53-11 1.73-13 

SSIDeCd) 

1.73-15 1.73-18 


3.13-12 2.93-15 

S3IDeC( 2) 

1 

b . 03-00 


1.93-16 
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¥g find that the computed order of S'5IDeC(l) is 8 
when m = 4, and 10 when m = 6. This is in accordance wita 
our theoretical result (3* 3.5). 

In practical situations the method of SBIDeC is very 
ef f icient. 

3 . 5 Oonclusion 

V'^ have studied the as37mntotic expansion for the 
global errors associatijd with the IDeC technique applied 
to quadrature methods for 'E'lBs. have shown that the 

order of convergence of an IDeC iteration is p'qual to th"" 
Order of the quadrature method on which it is based. This 
IS in accordance with the general result for IDeG proved 
by Stctter[l978]. 

Our study of IDoC has permitted the application of 
Richardson extrapolation on tho IDeC method and the presen- 
tation of tbe new acceleration t^-chnique SBIDeO together 
with its asymptotic expansion. A.s a consequ^nc'^, high 
order accuracy in the numerical solution of PIE (2.2.1) 
has b'^en obtained with only a modest amount of work. 

Another feature of SEIDeC is that its order of convergence 
IS unlimit'-^d s, wh-^reas that of IDeC is limited by the 
degr'^o of th polynomials used in the defect correction 
p rocedure. 



1/9 hav"* also described how SBIDcC is easily imple- 
as an automatic procedure. Kumerical example's 
d einonstrat'^ the efficioncy and th ? fast convergence of the 
new method. 

xln apnlication of SEIDeC is feasible only when the 
underlying IDeQ mcsthod admits an asymptotic expansion in 
powers of stepsize h. In the case of Volterra integral 
oq. nations 5 it is not easy to discover whether such an asymp- 
totic expansion exists for th" IDoC technique based on a 
given quadrature method. Moreover, our numerical results 
for thes= equations in Sections 2.3 and 2.4 do not exhibit 
any such asymptotic behaviour. Consequently, we have not 
considered SCIloC methods for Voltorra integral equations. 



CHAPTE R 4 


TK3 HITHODS OT ITERATED DEFERRED DEFECT CORRECTION 


. 1 I ntroduct ion 

In this chapt •’T vo pres'^nt a n-^w method of iterative 
iT'pro'^^ cment , -which is the result of applying the principle 
of Itora-*'ed Defect Corr-'ction ( iDoC ) on the method of 
Iterated Deferred Correction ( IDC ) , He call this ne-w 
method Iterated Deferred Defect Correction (iDDeC). The 
method IDD-^C is a generalization of the methods IDC and 
IDeC. It preserves tho good foatures of its constituent 
methods, and is free from their disadvantages. Our method, 
like its consul tu^^nts, finds solutions of increasingly high 
accuracy on th'^ same original grid (in contrast to th® 
Richardson extrapolation process which requires solutions 
also on finer grids). 

Section 4.1.1 contains a d-^scription of the method 
of iterated deferred correction as applied to FIE (2.2.1), 
and S'^ction 4. 2 a derivation of the asymptotic expansion 
of the method. 

In Section 4.3 ‘''re describe the method of IDDeC, and 
analyze its asymptotic behaviour. In Section 4.4 we tabu- 
lat ^ the numerical results obtained when IDDeC was applied 
to linear FIEs and to linear and nonlinear VIES, Our Diaplemen- 
tation of the IDDeC and IDC methods were based on the trapezoid 



64 


rule with Gregory corrocti)n and/or Sirrnson’s rule with 
corrr-c tion. 

4.1.1 Ite ra ted D-^f^rred Oorrection (l.DO) 

In order to explain the idea of deferred correction, 
use the linear Pr-^dholm integral eimtion (2.2.1) of the 
second kind: 

h 

(4.1.1) y(x) - jk(x,t) y(t) dt = f(x), a b. 

a 

¥g assume that the functions k(x,t) and f(x) possess conti- 
xiuous derivati^^^s of sufficiently high order, so that the 
uniq.ue solution y(x) exists and also has d'^rivativps of 
suf f ici ■'ntly high ord'^r. 

To introduce th-= discretization of (4.1.1), we let N 
be a positive integer, h = (b-a)/K, and x^ = t^ = a + ih, 

1 = 0(l)F. Purthirmorc, w"^ denote g( x^) by g^, where 
g(x) IS any function on [a,b], and repres'^nt the column 
vector (gQ, g^, ..., gj;r) g. 

¥? th^n consider th-^ Gregory's formula 
b 

(4.1.2) j'u(t) dt = h(i-Uj^ + u^ + ... + 

a 

+ h I Cg( 7^^ + (-Ij^ A®u^). 
s > 1 

with c^^ = -l/l2, ^2 = "•I/ 24 , = -19' 720 , c^ = -3/160, etc. 
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{ Bak-^r [1977] ) . We first neglect the second term on the 
right hand side of (4.1.2) and use the resulting trapezoid 
rule to write (4.1.1) as 

N 

(4.1.3) 7(x) - h B w k(x ,t ) y( t ) = f(x), a < X < h, 

tJ W ci 

where Y(x) is an approximation to y(x), and th^ w^ are the 
weights; w^ = % = t and w = 1 for 3 = l(l)F-l. The suh- 
stitution X = X|^, 1 = 0(1)N reduces (4.1.3) to a system of 
linear equations 

E 

(4.1.4) - h 4 k(x^,t^) = f(x^), 1 = 0 (l)H. 

0=0 

¥e denote th'^ solution of (4.1.4) hy 

now decide to include in (4,1.2) differences of 
order upto p (p< N) as Gregory correction to the trape-r^ 
zoid rule, and let 

. . P 

(4.1.5) = h Z OgL V®9^(T;b) + (-1)® ''So^(T;a)], 

s=l 

where 0^(Y;t) = k(x^, t)Y(t). Th^-n we solve the system 

E 

(4.1.6) Y^3_^^-h s w^ k(x^,t^)y^^^ = f(x.) + 

j=0 

and get und-^r certain conditions an improved solution 
This process of obtaining Y^^^ from Y^^°^ is called a 
deferred correction (see Baker [1977 ]) . 
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The procedure above may he repeated, and for 

— ( p ) 

r = 1, 2, improved so"^ ution v'^ctors T b-- found 

iterativply from 

(d.1.7) (I - = I + 

wh'^re K = [hw k( x , t )], and 

J ^ t] 

-( r ) 

Since th"' SGq.uencc of tho iterates ^ is not guarent^ed 

( n ) 

to converge, we take ^ ^ as the final approximation. This 
technique is now commonly called Iterated Deferred Correc- 
tion (IDG) (Baker [1977] ). 

1 variant of this method is an IDC in which equation 
(4.1.7) IS r^'^placed by 

(4.1.7^) (I - K;)Y^^^ = f + ). 

During the rth iteration of this method, the Gregory 
corr-^ction contains only diff-^rences of order upte r 
(a variable), and not upto p (a fixed constant) as in 
( 4.1.7). 

^ - 2 Asymptotic Bxpansion for IDG 

VJe shall now derive an asymptotic expansion for the 
global error in the solution produced by the IDC method. 

First, we derive, as in Section 3.T, an asymptotic 
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/ \ / \ 

Gxpansion for the global error y - \ where is the 

solution of 

( 4.2.1) (I - = ? 

and y is th*' solution of 

(4.2.2) (I - K)y = f + T 
in which a is givon by 

(4.2.3) T^S^^ly) 

b F 

= jk(x^,t) y(t) dt - h z w^k(x^,t^ 
a 3=0 

If wo choos--^ th'"' w corresponding to the trapezoid 

J 

rule, th^ quadrature errors howe the form 

(4.2.4) + A^(xj^)h^ + Ag(x^)h^ +' . . . . 

Now subtraction of ( 4.2.1) from (/|,2.2) yields 

(4.2.5) (I - E)(y - Y^°^) = P. 

Consoqueiitly, the expansion for 'cho global error is given 

ty 

(4.2.6) - y(x^) = + B^°\x^)h^ + ... 

Second, we derive a., error expansion for the solution 
obtained from after a single iteration of deferred 

1 2_ o 

correction given by 
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(4.2.7) (I - = f + 

whero is the Gr >gory corrsction to the trapezoid 

rule consisting of differences of order upto p, A rearrang-^- 
ment of the correction t-rrrs in (4^1*5) yields th^^ Lagrangian 
f 0 rm ”» 

p 

(4.2.8) O^P^(Y) = h L ^ [k( X , a+jh) Y(a+3h) 

J 

3=0 

+ k(x^j'b-;]h) Y(b-3h)]. 

The values of for p = l(l)4 are tabulated b?lov. 

J 


Table 4.2,1 Values of v^^p) 


3 

P 

0 

1 

2 

3 

4 

1 

- 1^12 

/ 1/12 




2 

- 1/8 

1/6 

-1/24 



3 

-109/720 

177/720 

-87/720 

19/720 


4 

-49/288 

77/240 

-7 '30 

73/720 

- 3 / 16 O 

Using ' 

( 4 . 2 , 6 ) and 

( 4 . 2 , 8 ), we 

can write 




(4.2.9) - C^^^^(y) 

= A*^2'^\x^)h^ + A\°\x^)h^ + ... . 
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where = !! 3+Dh.) B^*^^(a+ 3 h) 

;]=o 

+ k(x^, b-3h) B^°\b-3h)], n =. 2, 4, ... , 


If we now use th<^ Euler-Maclaurin series formula and write 
'i ■ as 

(4.2.4) X = S — ^ 

0=1 (23); 

where the B's are Be^'noulli constants ( Grant [1969 ]) , then 
the relation between the q.uadrature error and the Gregory 
correction may be expressed as in the following lemma. 


Lemma 4.2.1 

Ihe Gregory correction satisfies the eq.uation 

(4.2.10a) O^P\y) = T, R (x ) + 0(hP'^2), 

J- ' Pq 

where 


(4.2.10b) R (x^) 


Po+20 

(Po+2a-l) 

I [y_ 

3^0 (Po+23)' 


(PQ-tbjVl) 


J, 


With u(t) = k(xj_,t) y(t) and = p+2 or p+3 according as 
p IS eron or odd. (j-'bin followf- fro ' I’-^e of eourtionp 

(■^.1. f'^d (4o '’./!) 5 and formulae (nape 118, Baker 1877 ) for 
the derive liver in ter'^s of c^iff emcep of or(3 0 r p at inopt.) 

r/c observe th?t the 0 (b^^^) in (/, . 8 „irp) ig the error 
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In th ^ Is'ir'T.y if id giv=n by (4.2.4), the 

series B. ^^t) ^^7 wri+ten as 

Pq ^ 

( 4 . 2 . 10 c) R (x ) = A (x )h^° + A (x + ... . 

^0 ■‘- Pn 


But wG can write 


(4.2.11) - C^^P^(y) + C^^^^\y) 

= ( A^2° ^ + A'^^°\x^)h^ + ...) 

+ '^1 - R (x^) + 0(hP+2)^ 


from (4.2.9) and ( 4. 2.10a). This gives 
(4.2,11a) C^P^Y^'^b = -^ + 


where 

(4.2.11b) = {k^2°\x^)h^ + A^^°\x^)h^ + ...) 

- R (x ) + 0 (hP"^ 2 , 

Po 1 

Th'^refor'^ (4.2.7) can be written es 
(4.2.12) (I - = r + T + 

Subtracting (4.2.2) from tha above equation, we have 

(I - - y) = 44, 


and tho error expansion can b^ writoen in the form 
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(4.2 

.15) - y(x^) 

= 

+ . 




X ) -H 0(hV2), 



Po 

1 ’ 

wner 

stands for 

ohe same 

sines as R but with 


0 

P-i 

0 

a iff 

■’r^nt coefficients 






Third, give the nsyniptotic iD^hnviour r.n iterated 

A 

deferr''>d corrf^ction procodur' in the following theorem. 

? h'^ori^m 4.1 

L^t h-^ tne npproxunation to th'* solution of FIT 

( 4 . 1 . 1 ) "by th^ tr'^pezoid rule with stepsizo h, and let the 

_(v.) , , 

vectors , r = l(l)p satisfying 

(^. 2 , 14 ) (I - = I ^ ) 

ropr'^eent the successive approwiraa.tions hy the rr'-^thod of 
IDG using diff orences of 0 ’der upto p in the Gregory 
Gor'rection to the trapezoid rule. I^t p^ = p+2 or p+3 
according "S p is -^ven or odd, and suppose th^ p^th t-dari- 
vatives of the functions k(x^,t) y(t) ar= continuous in 
[e,b]. Then the asymptotic expansion for the global error 
IS given by 

(/u2.15) ty’ - y(x^) = ^ ^ __j 

( IT ) 

+ S (x ) + 

Po 

where ) is of th^ form R of (4.2.10). 

Po ^ Po 



P roof 


Vo hav- already shov i th'' t xii? thoorera xs truo for 
r = Supposo it IS true for r = 1 , 2 , rQ(r^'<p). 

ih^x- (<, 2 , 15 ) IS true, in p^^rticular for r = r^. Th'^t is. 


.( ) 
1 




r +2 ( r ) 


( 4,2 16 ) ° - y( X ) = (B,, ° (x )h ° ° (x, )h ° 


r +4 


4 


+ S ° (x ) + 


Me shall now proTo that ( 4.2.15) is tru'^ for r = r +1. 


Th^ oq.ua tion (4.2,16) i^npli-'S 


^ ) - C^P\|) =(^2"° (x^t 


)ii " (xjh 


^ 0+5 + 


+ ' 2 ^ ° (x ) + 0 (hP'^'^), whore 

Pr, ^ 


A 


(rJ 


n 




( 


[k(x^,a+;^h) (a+D^) 


3=0 


/ 


+ k(x^,h-Oh) ° .'h- 3 h)j 


^ 4 ) 

and T (x ) has the s'^rao fori^. as R 


p 


1 


Po 


Again using (4.2,10) wo can writ*' 

_ Jv-i) 

(l-k)(Y -y) = 0 , whore 

(r + 1 ) (r ) r +3 (r ) r +5 

= (A2 ° (x^)h ° + ^ (x^h 4- , ^ 

+ ) - R (x ) + OlhP-^^). 

Pn ^ Pn ^ 
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' 3 Itorntod Deferred Defec~ Corrocticn (IDDgG) 

. 3. 1 D'^ scri-ption of the m'3t/iod of IDI 

'./e now doscrilos an application :f the method of IDeC, 
which wn,s ?xplained in Section 2.2, o.i the mf^thod of IDG. 

■'•■'c assume that ¥ is a multiple ')f m, where m is <=v^n 
and greater than max[p, 2., with the meaning of p as in 
5). L"^t denote the IDO solution of the 

FIB (/i.l.l). 

'Jo dD.Tide [a,h] into subinter^'als I^ = [x( ^4.1)1-15 L 

3 = l(l)lT/m. For each ;], lot th» polynomial of 

degree m interpolating the m-f-1 pofits 
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( 4 1 3K1. We then define, for x in [n,!:], the 

interpolating function H x) as the collection of poly- 
nomials such that 

J r, m 


(4.3.1) 


Q 

'in 


(x) = 


P*^^\x) 

3,1^ 


for X e I 


3 


= l(l)N/m. 


Consider now th^- problem 
b 

(^i«3.2) y(x) - Jk(xs,t) y(t) dt = g^°\x), 

a 


a < X < b 


where 


b 

g^°^(x) =P^°^(x) -/k(x,t) P^°^(x) dt. 

a 

Since IS a good approximation to y(x), the 

expression - f(x) is ■expected to b'^ small and thus 

(4.5.2) IS a n-^ighbour ing probloia (KP) of our original 
problem (Op) (/j-.l.l). The -^xact solution of this new 
problem is th-^ known function l’^°^(3^). solv-^ this UP 

using the sejae m'^thod by which W'^ solved our Op, and denote 
the resulting solution by I7 :s known error 

r^°\x^) - baken as an estimate for ohe unknown 

orror y(x^) - in the IDG solution of our OP. 

Thus, to th= true solution y(x), th= approximation Z^ 
given by 


2U,p) ^ 2U,p) ^ 


( -I.Z, 4) 
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IS nor 3 accurate than 

rep oat the process above, and find improved 
appro-xm tions s = 0, 1, 2, giv^n by 

(^h^.5) zU+l-p) ^ z(^o,p) +p^s)(^^) _ 

whc-rr' the function def^n^d analogous to the func- 

tion P^^*^ x) of (4.3.2), int ■'rpoD.ates the points 

( und is th^' IDC solution of the following 

np; 

b 

c 3 . 5 ) y( x) - j k( X, t ) y( t ) dt = g^ ® ^ ( x) , a < x < b, 
a 

b 

with g^"^(x) =P^®Hx) -/k(x,t) dt. 

a 

Thus, storting with the solution produced 

by Lh TPC^ we find its iterative improvements by the 
metbod of IDeC, each iteration of which involves an IDC i 
Hence we call the '^ntiro t'‘chniq.u'' the method of Iterated 
Deferred D^f -^c c Correction (iDDeC). 

» 3 • 2 j-Ls yrTot o tic '^xiQ' ' ^nsion for IDDeO 

1' 3 r 'Call the st-^ps carried out during th'' first 
^tcra'^ion of our IDDeC. first, the solutions of our OP and 
HP are obtained as depict-^d ir Tabl^ 4 . 3 . 1 . S-=cond, the 
improved approxirnation of is computed 



using ( 4. 3. 4) . 


“t ^ 

VI 0 now proco'^d to dPx’iv3 an 

i7 * 

che global orror in tb.'^ first it^rat-^ •'-'-I 
Theorom 4,1, the solution foun^ ‘ 

ScatJsfios th : op nation 

( 4 . 3 . 7 ) - y(x^) == • 

Po ^ 

1 * ** 

whoro S^p\x ) = )b^° + 

i^o Pq '^o 

is either -th-^ intogf^rs p+2 or p+3 wb i * 

Tn blv. 4.3 . 1 iLey st^ps_iJ.-- 


Probl'^ia Method. dystem solved 

cr-^pozoid = f 

Op 

IDO with (I-K)Y*^^^ = T + C^P^T' 
r=l(l)p 


KP 


trapozoid (l-K)y^°’°^ = g 


IDG with 


r=l( l)p 
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7 3 now nssunip p is '-von. Since int''rpolat ^s 

m 

the points ( Ti 3 h-nve 

(4.3 8) pj^°\x) = jix) + + ... 

+ 4- 0(hP^-2), 

provided p+2 < n+2, nnd th.'- ■^rror duo to interpolation is 
"'ssuirod to bo included in the 0(hP'*'^) t^rm. 'Th.or^for'? th'^ 
ppror* 


in th - solution of NP ( 'I-.;.-?) can bo writt'^n as 
(4,3 9) T^°' = + 0(h5+t. 


Analogous to th^ oqtations (4.2.5) and (4.2.6) in tho 
case of Op, tho equations for our Np are 


( 4 . 3 . 10 ) 


(I . 


ind 


(4.3.11) 


..( 0 , 0 ) _ ^ + 


+ B'°)(xphP+2 - B^^(xphP+'^ + . 


Therof or^j 
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(4.3,12) » l[°hx^)h^ 4 ... 

+ A^“)(xj_)h5+3 + D'(«)(:Ci)hP+'‘ + ... . 

But we hav-' 

(4.3.13) + 0(hI’+2) 

and ( 4 . 3 . 12 ) b'^coiaes 

(4.3.14) = (A^2®\x^)h^ + ,,. + A^°^(x^)hP‘^^ 
+ D-(^(K,)hP -5 +...)+ 4=) - 

An application of IDC to FP rise to 

(4.3.15) (I-K)y^°’^^ = 

where 

= (A^°^(x^)h^ + ... + •^p°2( **■ T ' 

+ D'^^{xphP+5 + ...) + S^_^2(x^) + 0(hP+2). 
Moreover, 

Ci-k)(y^°-i) - 5^°)) = 

and 

(4.3.16) - P^°^(x^) = (B^2^^(Xj_)h^ + ... 

+ B(i)(xphP-5 + +•••)+ 
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After r applications of deferred corr<=ction, we have 
f 4.3.17) + ... 

+ B^j)(xphP+2+^ + D^P)(xpilP+‘^+^ + ...) 

+ S^P)(x^) + 0(hP+2)_ 
provided p+/l-+r < iri+2. 

Vg note that tho torn' OCh^"^^) in (4.3.17) is th^ 
as th'^t in (4.2.15), except for th''^ error due to interpo- 
lation. Using Zadunaisky's principle, we therefore h'^ve 
bhe result; 

if 2p+4 < ni+2, then 

(4.3.18) = y(x ) + )h^P^^ 

1 1 P -T 1 

+ + ... 

P+3 1 

A si'Milar result holds when p is odd, and the follo- 
wing theorem is true, by induction. 

Theoreiii 4.2 

Under the conditions stated in Theorem 4.1, for each 
s = 0, 1, 2, ... the solution produced at th^ end 

of the sth iteration of the IDDeC, the UPs of which involve 
polynoinials of degree m, has an expansion of the 


f orm ; 
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(4.3.19) Z 


( s.p) 




+ B 


(s,p) 


( 


+- 


b(s,p)( 

\+2 ^ 


x^)h' 


g+2 


+ 


■where 


(4.3.19a) q = miii[( s+l)(p+2) , ra+2]. 

Rea^rks 

Fro. I Thooren 4.2 above, it is se^n that '''/’hen s=0 our 
IDDeG reduces to an IDG, and the aaxin'uin order of accuracy 
obtainable is p+2 (cf. Theorem 4.22 in Baker [1977 ]) . On 
the oth-^r ha,nd when p=0. our IDDeC degenerates into an IDeC, 
end th"' piaxiaua order of accuracy q attainable after 
s iterations is niin[2(s+l), a+2]. In general, if the ordor 
of the underlying method is n, th ^n q = min[(s+l)n, m+l]. 
Thus, our r 'suits agreo with those concerning the iDpO 
othods applied to differential end integral equations, given 
respectively by Frank [1975, 1976, 1977] and Sathiaraj and 
Sanka,r [198^] (se? Theoren 3.2.1), and IDDeO is indefed a 
griioralization of both IDG and IDeQ. 


4.3.3 Imnlementation of IDDeO 

Vo givf’ below an outline of the IDDeC as applied to 

( ^ 1 . 1 . 1 ). 


Algor it hi: IDDeC 

1. Input a, b, IT, h, p, and m. 


2. Find the solution of the Op by the use of a basic 
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discretization rule nnd then find the solution Z. 

1 

hy IDG using upto the pth order differences as 

Op 

correction teri'-is to the rule. Copy in . 


For s = 1, 2, . . , , do 

5.1 Find 3 pipcewis^ int-^rpolating polyno’^i‘^1 
P^(x) of degree m fitting the data points 
(x , Z ). 

3.2 Sol'V'-^ the FTP first by the basic rule and then 

by the sane IDG method used in Step Let 


its solution be Z‘ 


FTP 


3.3 Improve the solution, using 

Z. 


Z?^ + (Z. - Z^"^) 
1 1 1 


4 , Output Z , 1 = 0 ( 1 )N, 


In Step 3.2, the solution of NP involves the calcu- 
lation for 1 = 0(l)]I of the values 


b SM d 

/,3.20) I = jk(x ,t)p '' t)dt = L jk(x ,t)P (t)dt. 

I-I -L Hi ^ J 

3=1 a 


For each 3, we find the polyncnial P a.nd then 

Hi ^ J 

c 1,010 tho 001 responding values of 

b 

/k(x ,t)P (t)dt 

-L Li , J 

a 

in on array called YALUTT (1,3), 1 = 0 (l)U, 3 = l(l)l\T/m„ 
Then the values are found using 
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N/m 

= £ VAIUTT 
3=1 

¥e mado use of the HAG- routines while calculating the inte- 
grals in ( 4.3.20). 

4 . 4 Ap-plications to Integral Equations 

This section deals with applications of the methods 
of iterated deferred correction and iterated deferred defect 
correction to il lustrative examples in linear Fredholm, and 
linear and nonlinear Volterra integral eq. nations, 

4.4.1 Linear Fredholm Equations 

¥e implemented the method of IDLeG as applied to 
(4.1.1). The techniques of IDDeC and IDG were based on 
quadrature methods which make use of the trapezoid and 
Simpson’s rule. 

(a) The IDG methods based on the trapezoid rule 

We give results from th'^ computational tests performed 
using 16-18 significant digits. Problems 1 to 4 of 
Section 2.2,2 were solved by IDCs based on the trape/.oid 
rule teg.-^ther with th^ Gregory correction and by IDDeGs 
generated by the above IDGs. 

¥e denote by II)C(p,s) the IDG method involving s 
iterations of deferred correction with upto pth order 
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correction terms to tii? trapezoid rule. Similarly the 
notation IDDr'C(p,s) stands for th= IDDeC method which 
involves s iterations of defect correction end which is 
based on the method IDC(p,p). 

fhe problems wore solved by those methods for different 
values of p, s and N; th'^ de;STee m of the polynomials used 
in IDDeC was chosen to be 4, 6, and 14. In the 

following tables we display the maximum absolute errors in 
the computed solutions and the errors with asterisks (*) 
indicate that they remain the same for the subsoq.u'’nt 
iterations/ corrections. 

V/g tested the performance of the two types of IDC 
defined by (4.1.7) and (4.1.7*) and found them essentially 
the sama (se'' Th'^orem 4.22 in Baker [1977] ) . Tablo 4.4.1 
displays the results for Problem 1 by these two types of 
IDO with Gr-^gory correction where q. is a constant p 

or a variable r, and those by the corresponding IDDeC 
m'-^thods. ¥e wish to add that the methods IDDeC(4,l) and 
IDDeC(3,l) yielded the same results as IDDeC(2,2), whether 
q IS a constant oravariable. 

We also tested an IDO method in which th'^ first 
deferred correction involved differences of first order 
only, but th-^ second correction involved differences of 
upto the third order. The IDDeC method based on the above 
IDC was also implemented. The results of these methods are 
listed in Table 4.4.2, 
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Table 4.4.2 IDC, and IDDeC with different values of m 



N 


IDC 


f 

t 

1 

t 

IDDeC 

(1,0) 

(1,1) 

(3, 2) 

1 ^ 

.1 

Iteratron 1 


12 

3 .9E-3 

3.5E-4 

1.3E-6 

1 

! 4 

1.2E-8* 






! 6 

1 

1 

5.5E-11* 

Problem 1 

24 

9 . 6E-4 

4.5E-5 

4.1E-8 

1 

1 

I 4 

1 .9E-10* 






I 6 

2.1E-13* 






: 8 

1 

2.7E-16* 






1 

1 _ 



48 

2.4E-4 

5.7E-6 

1.3E-9 

1 

1 

I 4 

2.9E-12* 






1 6 

8.4E-16* 






I 8 
\ 

1 

1 

4.3E-18* 


12 

2,7E-3 

3.0E-5 

3.5E-6 

i'' ' " 

\ 

1 

I 4 

3.4E-7* 






! 6 

1 

1 

1 

1.2E-8* 

Problem 2 

24 

6.8E-4 

1.9E-6 

5 . 6E-8 

1 

1 

1 

1 4 

5.0E-9* 






1 6 

4.5E-11* 






! 8 
i 
) 

4.6E-13^ 


48 

1.7E-4 

1.2E-7 

8.8E-10 

pM— — . 

1 

1 

1 4 

8.0E-11* 






1 6 

1.7E-13* 






1 8 

1 

1 

1 

» 

.j 

4.4E-16* 
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Table 4.4.3 IDO with q z= r = 1 2, and IDDoO with m = 4 


Problen F 

IDC(1,1) 

IDDeG( 1,1) 

1 

IDDeC (1, 2) 

IDC( 2,2) 

IDDeC( 2,1) 

12 

3.5D^4 

3.8E-8 

1.2E-8'^ 

1.2E-5 

1. 2E-8 " 

1 24 

4. 5B-5 

6.23-10 

1.9E-10* 

1.2B-7 

1. 9S- 10* 

48 

5,7B-6 

9 . 8E-12 

2.9E-12^ 

5. 6E-8 

2.9E-12 ^ 

12 

3. OB-5 

3.3B-7 

3 . 4 B- 7 '' 

5.8B-5 

3.4E-7* 

2 24 

1.9B-6 

5. IE-/ 


: 3.7E-6 

5.1E-9* 

48 

1.2B-7 

7.9B-11 

8 .OE- 1 / 

2.3E-7 

S 

8 .OE-II"' 

Fext , 

we list 

in Table 4. 

4. 3 th'' results by th' 

IDO 

and IDDeC nii 

ethods for which the 

maximun ord 

■^r of diff 

"■rences 

q in the G-r 

"■"gory correction term is taken to Toe 

th'’ 

ij/(tera^tion 

nunber, 

r(r = 1, 2) 

. Wo find that even for the 

snail veiu-^ 

4 of rn the IDDeC produces v^ry 

good r'-'sults. 

The nuiiorical 

results for 

the two typos of IDC 

methods, 

with q = p 

or ry are 

essentially th-^ same. 

This IS true 


also for th'^ corresponding IDDcC nethods which are generated 
by thi^so IDOs. Th-^refore we list in the succeeding 
Tables 4.4.4 to 4.4.8 the results only cf the first versioa 
of th ^ IDO rn'‘thod, end those of the corresponding IDDeC 


mc'thod. 


Table 4.4,4a IDC and IDDeC methods wxth p 


1 


Problem 

m 

N 

IDC 

IDDeC 

1 

Iteration 

2 

number 

1 

6 

6 

2.7i-3 

3»lE-6 

1.7E-8 

1.4E-3 



12 

3.5E-4 

5.0E~8 

6.2E-11 

5.5E-11 



24 

4.5E-5 

8. lE-10 

2.3E-11 

2.1E~13 


8 

30 

1 . 2E-3 

5.6E-7 

2, 3S-10 

1»5E-11 



16 

l,5£-4 

9.1E-9 

5 . 13 

1.5E-14 



32 

1.9E-5 

1.4E-10 

1,1E-15 

1.5E-17 


10 

10 

6.1D-4 

1, 5E-7 

3 ,6E-1] 

3 .9E-15 



20 

7.7E-5 

2.4S-9 

7.4E-14 

5.3E-18 



40 

9. 8E~6 

3.8E-11 

1 .5E-16 

4 »3E-18 

2 

6 

6 

4, 8E-4 

4 . 8E- 6 

4,3E-6 

4 .5E-6 



12 

3 .OE-5 

1.7E-8 

1.2E-8’' 




24 

1.9E-6 

7.8E-11 

4.5E-11=^ 



CO 

8 

1.5E-4 

3„3E-8 

3 .9E-8 




16 

9 . 4E-5 

5,4E-10 

2 .8E-11 

NC 



32 

5.9E-7 

4.4E-12 

2.6E-14 


3 

6 

6 

2.4E-3 

5.7E-5 

1 . 2E-6 

2.3E-8 



12 

3.3E-4 

1.2E-6 

3 .5E-9 

l.lE-11 



24 

4.4E-5 

2.2E-8 

8.9Ew12 

3.7E-15 


8 

8 

l.]E-3 

1.2E-5 

l.lE-7 




16 

1.4E-4 

2.3E-7 

3 .OE-IO 

NC 



32 

1.9E-5 

4.0E-9 

7.2E-13 


4 

6 

6 

3.0E-3 

8.8E-5 

1.6E-5 

1.8E-5 



12 

3.7E-4 

1.4J1.-6 

7.3E-8 

7.7E-8 



24 

4*7E-5 

2.5E-8 

2.9E-10 

3 .OE-10 


8 

8 

1.2E-3 

1.4E-5 

2.9E-7 




16 

1. 6E-4 

2.6E-7 

5.4E-10 

NC 



32 

2. OE-5 

4 . 6E-9 

l.OE-12 



>ote: NC stands for 'not computed’ 



Table 4.4 » 4b Computed orders for the iterates of IDDeC 
with p = 1 


Problem 



IDDeC iteration 

m 

0 

1 

2 

0^ 


6 

3 

6 

8 

8 

1 

8 

3 

6 

9 

10 


10 

3 

6 

9 

L 


6 

4 

6+ 

8 

8 

2 







8 

4 

6+ 

10 

1>1C 


6 3 6 8+ 1 0+ 

3 

8 3 6 9 HC 


4 


6 3 6 8 8 

8 3 6 9 NC 







Table 4.4.4b Computed o rders for the iterates of IDDeC 


with p = 1 






IDDeC 

Iteration 

Problem 

m 

0 


1 

2 

3 


6 

3 


6 

8 

8 

1 

8 

3 


6 

9 

10 


10 

3 


6 

9 

L 


6 

4 


6+ 

8 

8 

2 

8 

4 


6+ 

10 

l^C 


6 

3 


6 

8+ 

10+ 

3 

8 

3 


6 

9 

NC 


6 

O 

O 


6 

8 

8 

4 

8 

3 


6 

9 

NC 
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Tablos ,4.4.4 to 4.4.6 give the n®-.xiivuiia errors of th"^ 
methods IDC(p,p) end IDDeC(p,s) for various velues of p, 
s ond K when applied to Prohlems 1-4. 

Table 4.4.4s gives th-^ '^rrors in th-’ solution for 
Probl^’-ns 1 to 4 by th^^ IDO a'^thod with p = 1 and by the 
IDDeC method with differ'-^nt values of m, and Table 4.4.4b 
depicts th''^ computed orders of convergence for the iterates 
of tho ILDeC methods. These computed orders are in accor- 
dance with the theoretical orders predicted by tho 
eq.uation (4.3.19a). In the table a number with n '+' sign 
moons the computed order is more than that number. 

Tabl' 4.4.5 portrays th^ validity of order relation 
(4.3.19a), for m = 6(2)10 and p = 2(l)4 in th'’ case of 
Problem 1, and Table 4.4.6 for n = 6 '^nd 8, and p = 2 in the 
cas'’ of Problems 2 to 4. 

Tables 4.4.7 and 4.4.8 display the results by the IDO 
and IDDoC m.^'thods for diff-^rent value's of p, and reveals 
th"' power of IDDeC over a wide rang^' of m and also for 
relatively small values of K. 

Our computational t'^sts show that the IDDeC method is 
sup'^rior to tho IDeC a,s regards accura,cy and cost. 



Table 4.4.5 The errors for Problem 1 by IDC(p/p) and IDDeC(p/l) 


where p = 2, 3^ 4 . 



N 

m = 

t 

1 

■ 

■ 

8 

1 

1 

I 

m - 

10 



IDC 

IDLeC 

T" 

I 

IDC 


IDDeC 

1 

1 

I 

1 

IDC 

IDDeC 

p=2 

m 

2S-4 

lE-8* 


6E- 

5 

lE-11* 

1 

1 

1 

1 

1 

1 

1 

3E-5 

4E-13^^ 


2m 

IE- 5 

5E-11* 


4E- 

6 

1E-14>- 

t 

1 

1 

t 

2E-6 

4E-16* 


4m 

9E-7 

2E-13* 

-i- 

3E- 

7 

lE-17-^ 

» 

1 

t 

1 

1 

lE-7 

4E-18* 

Order 


4 

8 

M 

4 


10 

"HT" 

I 

1 

1 

i 

4 

8 

p=3 

m 

6E-5 

lE-8* 


IE- 

5 

2E-11* 

1 

1 

1 

1 

f 

5E-6 

lE-14’^ 


2m 

2E-6 

6E-11* 


4E- 

7 

2E-14-* 

1 

1 

f 

1 

2E-7 

4E-18* 


3m 

6E-8 

2E-13* 

\ 

-4- 

lE- 

8 

2E-17’- 

1 

1 

1 

f 

1 

5E-9 

4E-18* 

Order 

5 

8 

«JL 

5 


10 

1 , 

1 

f 

f 

1 

5 

10+ 

p=4 

m 

2E-6 

lE-8* 


5E- 

■ 7 

2E-11* 

'"’i 

1 

1 

1 

t 

2E-7 

lE-14* 


2m 

6E-8 

6E-11* 


lE- 

■8 

2E-14*- 

1 

: 

1 

f 

3E-9 

3E-18* 


3m 

lE-9 

3E-13* 

j 

2E-10 

2E-17* 

1 

1 

f 

1 

1 

5E-11 

4E-18* 

Order 

6 

8 

-1 

6 


10 

j 

1 

i 

I 

1 

6 
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Table 4.4.7 IDC and IDDeC methods with p = 2 and 4 


} 

— _ 1 


m = 

12 i 

m = 

14 



P 

N ; 

IDC 

IDDeC ; 

IDDeC 




- - -1 

(p/p) 

(p,l) 

(p/2) 

(p.l) 

(p,2) 

Problem 

2 

2 

m 1 

5.8E-5 

1.8E-12 

1.5E-12* 

5.9E-14 

6.3E-1' 




2m 

L 

3.7E-6 

3.7E-16 

6.6E-17* 

6. 5E-17 

3 .9E-1! 

1 

Order [ 


8+ 

14+ 

8+ 

L 



4 

m 

2.0E-6 

1.5E-12-^ 


6.3E-15-^ 





2m 

3 .3E-8 

6.8E-17* 

. . 

1.5E-1S 

3.lE~lt 

1 

Order I 

-- - - 


12+ 

14+ 

12+ 

L 

Problem 

3 

2 

m 

1.4E-5 

1.9E-8 

1.4E-11 

6.0E-9 

2.4E-1J 




2m 

8.8E-7 

9.6E-11 

5.0E-15 

- . - 

2-9E-11 

8,3E-l^ 

Order 




8 

12 

8 

12 



4 

m 

5.3E-8 

1.8E-11 

7.2E-16 

3 a5E-12 

7.2E~11 




2m 

l.OE-9 

8.7E-15 

3.9E-18 

1.5E-15 

7.8E-ie 

Order 


12 

L 

12 

L 

P roblem 

4 

2 

m 

1.5E-4 

4.9E-8 

3.5E-11 

1.5E-8 

7,3E-12 




2m 

l.lE-5 

2.5E-10 

1»5E-14 

7.7E-11 

2.6E-15 

Order 


8 

12 

8 

12 



4 

m 

4.9E-7 

3.2E-10 

6. 2E-12 

5.6E-11 

2,8E-14 




2m 

1.2E-8 

l.lE-13 

5.2E-16 

1*9E-14 

4.8E-18 

Order 


12 

14 

12 

L 
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Table 4.4.8 Errors in so lutions by IDDeC 




m 

i 

4 : 

m = 

12 



IDDeC 

IDDeC 

Probleni 

N 

(2, 1) 

(2, 2) 

. ... 

(3,1) 

(3, 2) 


m 

S.4E-6’' 


2.3E-17 

8.7E-18 

1 

?m 

1.4S-7* 


2428-18 

4.3E-18 


m 

3.4S-4* 


1.5E-12 

1.5E-12 

2 

2m 

4.0D-6* 


6.9E-17 

6.7E-17 


m 

5.6E-5 

1 4 8ij - 6 

8.6E-10 

1.5E-13 

3 

2m 

3 , 8E- 7 

1.2E-9 

1.2E-12 

l.lE-17 


TV 

1.5E-3 

1.4E-3 

3 . 5E-9 

5 4 5E~12 

4 

2m 

2.4E-5 

2.3E-5 

5.2E-12 

4.8E-16 


(b) The IDO methods based on Simpson" s formula 

Just as Gregory’s forniula represents the trapezoid 
rule plus a correction, the following eq.uation represents 
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Sanpeoia's rulo with i corr^'^ctlon (so'-'' Fox[l962]): 

( . 4 .> 1 ) ' udx = •^ ( + 4u^ + 2 U .2 + . . . + ]_ ^ ^ 

^0 

N-1 

_ l3-„ V a4 . 

90 ^ + .oc, 

0 =1(2) ^ 

IM IS iv-'ii and 

(‘^.4.2) 6 '^u^=u,^-4u -|4-6 u -4u -,+u 

0 0^2 ;]hl 0 0~1 0-2 

L^t us apply th'^ quadraturo mtliod h^^sed on Siupson's 
rul ' to th FIE ( 4 . 1 . 1 ), and obtain tho solution ^ fro.i 

(4.4.3) (I-K)Y^°^ = f 


TTh.-'r^’ L is an-ilogous to th > u'^trix in (4,1.7). 

'.nth bhv-^ notation of S '■c tion 4,1, a sin^"; 3.pplica ' lOii 
of doforrod corr'^ction gives froi^ 

(/,u4,4a) (I~K)Y^^^ = f + wh-rp 


4.4.4bj = 


90 


[6^0^(Y; t^) + 6^0^(Yj t^) H- 




w_'uh tj =. k(x,t)Y(t) . 


"to nox-" consid-^r an IDG nethod d'^fin^d by 


(4,4.5) (T-K)Y' 


( r) 


+ U 


( 4)^y^ ) 


r = 1, 2, 



96 


a knowledg'^ of k(x^yt)Y(t) and f(x) for points t = x = x_^ 
and which, are outside the basic interval. But then 

we can comput-^ th-^ values outside the range with the use 
of (4.1.1) or 

(4c4.7) T^^\x) - S(x) = f(x), 

where S(x) represents the Simpson formula: 

S(x) = |[k( x,Xq)y[j°^ + 4k(x,x^)Y^°^ + 2k( x, X2)Y^2° ^ + ••• 

+ 4k( X, ° 1 + X, Xjj ) Y^^ ° ^ ] 

Then the use of (4.4.4a) gives Y^^\ 1 = 0(l)Er. 

Again th^ corr-^ct-^d values and niay be foxmd from 

(4.4.8 ) Y^^Hx^) - S(x^) = f(x^) + 

with 1 = -1 ana H+1. Similarly the corrected values 
1 = -l( l)N+l may th'^n !<= found using (4.4.5) and (4.4.8). 

Our IDDeC proce-^ds by taking =Y^^^ and 

finding the solution of th^ sth hP as we have don? in the 
case of Op. As for the ^valuation of g^ ( x) at x = 
and we assume \ ^ = -1? N+1. 

llPISYJL4P_ 3-1 r-^sults 

Problems 1 to 4 of Section 2.2.2 were solved by the 
IDG and IDDeC methods based on Simpson s rule with 
correction (4.4.1). In Table 4.4.9 wo list the errors in 
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In our implementation of IlDeC, we solved the NPs by 

( s ) 

taking ^(x^) = f(x^), i = -1, N+1. Mor ^over, we assumed 

= y(x^)^ 1 ^ N+1. 

Our computational results indicate that the order of 
the IDC(2) method is 6 as expected. The IDDeC method 
produces highly accurat'^ r^^-sults but the order of accuracy 
IS not as high as predicted by the th.-^oretical analysis. 

JO note that the IDDeC method based on the trapezoid 
rule with coirection upto fourth order differences as well 
as the on’^ based on Simpson’s rule with a correction (4.4.1) 
IS expected to produce solutions of the same order, namely 
12, 18, Or m+2. However, we would rath-^r choose th^ 

former becaus'^ of the simplicity of th^ trapezoid ruli^ and 
b'caus'^ of th"' implementation difficulties associated with 
the latter IDDeC. (Compare results for m = 12 in Tables 
^^4.7 - 4.4.9. ) 

4.4.2 Li nee.r Volterr a Integral Equ ations 

In S-^ction 2.4 w'' consid'^^ red the quadrature methods 
for linear Volterra integral equations of the second kind 

X 

y( x) -/ k( X, t) y( t) dt = f ( x) , 0 < x < a. 

0 


(4.4.9) 
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Fov wr- shall descrihs two approachos of itaratad def'^rr^d 
correction. (,ID0) hased on a Quadrature method. ¥e tested 
tho p'^rf orriances of th-^S'^ IDCs and the gon'^’rated IDDaCs. 

J'FlLf * ^iL ^ oj^3 ?r»ctiorL - local aunJ^oach 

A process of def 3rrpd correction based on Baker[l977] 

which vr call a 'local' approach is now applied to ( 4 , 4 , 9 ). 

rh 

For th'' integral j'u(t)dt, Gregory's formula may be written 

0 

a s 

rh p 

(/|.^„lo) j u(t)dt = T + h Z c [ v'^n(rh) + (-1)® '^®u(0)], 

0 s=l 

wheio T represents the trapezoid rule approximation with 
stepsize h, and the other expression on the right hand 
siQp is th-2 Gregory correction analogous to (4,1.''). 

TTst ig th-^ above eQuaiion (/1.4.IO), we write the 
c Qua lions 

Y(0) = f(0) 

--|hk(h,0)T(0) + [1 - ihKh,h)]Y(h) =:f(h) + C3_ 

( 4 , 4.11 ) 

-■^hk(2h,0)Y(0) - hk(2h,h)Y(h) 

+ [1 - |hk(2h,2h)]Y(2h) = f(2h) + C 2 

whore C^, r = l(l)F r'^prosents Gregory cor'^ection in 
(^,4.10) with u(t) = k(rh,t)Y(t'. 
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In 1113 rth slagr of th'^ local approach, of dof^^rred 
correction, we calculate T{ rh) in "terms of pr'^"viously 
accept-d approximations to Y(0), T(l), YCF^lh). The 

rth stage is initiated hy computing Y^°\rh) from the 
following formula in which is taken as zero. 


r-1 


(4. ^.12) [1 _ ■|hk(rh,rh) ]Y^°\rh) = f(rh) + h £ k(rh, oh)Y( jh) 

D=1 

+ ■^hk(rh,0)Y(0) + 0^°^ 


The finite differences of the values k( rh, 3 h)Y( jh) , 

3 = 0(1) r-1, end k(rh, rh)Y^°^ can be used in the Gregory 
corr'"'ction to th'^ trapezoid rule in order to obtain a 
correction term 0^*^^. th^n find a corrected value 
Y^^\rh) satisfying 

(4.4.13) Y^^\rh) = Y^'^^rh) + C^°V[3. - IhkC rh, rh) ]. 

The process can be repeated, using Y^^^(rh) to 
compute and th-^reby Y^^\rh), and so on. The process 

mey b? continued for a fixed number s of iterations or 
until Y^®\rh), for soni'=' s, agrees to the required accuracy 
with Y^^'‘^^(rh). When Y^®^(rh) is accepted it is taken 
n s Y( rh) in th-^ formula used for calculating the subsequent 
values. Th-' above m'^thod of IDG involves a local correction 
since Y(rh) is corrected before Y[(r+l)h] is. 
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P eforrod O o rraction - Glolal A-p-proach. 

W'e now introduce a ^lolal approach, of deferred 
correction. In this approach we first find th® solution 
Y(^h) = Y^°\rh), r =0(1)15 by the quadrature method (4.4.11) 
and th'ii compute all the correction terms r = l( l)K 

using" th'^ expression for Gregory correction. Then we 
cornputo 9 13. fch^ corrected values Y^^\rh), r = l( l)]S from 

(4.4.13). 

Using T^^^(rh) wo compute and thereby Y^^Hrh), 

and the process is repeat'^d adopting a stopping criterion 
similar to th'^ one employed in the local IDG. 

This approach is called 'global', b^^cause correction 
is made to all values of Y simultaneously. 

N um er real re sults 

The two approaches of IDG described above have been 
ir'p] -m'^nt-'d together with the iterated deferred defect 
co-'rection 'a'"’thods based on these approaches. The test 
runs hev'"^ been performed in single precision ( 9— H decimals) 
on oh.'" IGL 190 /jS computer at BKU, Varanasi. ¥'=■ consider 
'ipair til"' l-'ct problem given in Section 2.4. More test 
r'-fn’i Is Mr"' included in Chapter 8. 

1 li tli^ t'^st runs, we took thv"’ number of correction 
terms p Oi equation (4.4.10) as 1 to 4, and the number of 
Iterations s as 1, 2 or 3. The basic interval chosen was 



-4.10 P esu 1 ts for local IDC and IDDeC 


’\rr, 

IDC 

1 

IDDeC Iteration 

2 3 

number 

4 

5 

i- 

1. lE-1 

3,3E-.3 

5.9E-4 

7,9E-4’‘ 




3 .7E-2 

7.95-4 

5,lE-5 

3.9E-5 

2.1E-7 

1.9E-8 

L? 

3 

2.51i-4 

1.6E-5 

l,0E-6 

6.4E-8 

4.3E-9 

1C 

2. 3E-3 

l.OE-4 

3 .7E-6 

3.8E-6 

9.2E-6 

5.0E-7 

/ 

1 

3 . CL- 2 

1.3E 3 

b.3E-4 

5.9E-4* 



e ) 
o 

. IE -3 

1.2E-4 

2,0E-5 

3.6E-6 

6.2E-7 

1.2E-7 

12 

1. 2E-'. 

3 .5E-5 

9,7E- S 

2.7E-6 

7.8E-7 

2.2E-7 

16 


1. 3E-5 

6.8E-5 

4.0E-6 

2.1E-6 

8.3E-7 

i 

2.9L-2 

5„3E-4 

3 .OE-4* 




f; 

4.7L-3 

2.3E-4 

2-0E-5 

1.8E-6 

1.5E-7 

l.lE-8 

1? 

1 , 3E-3 

1. lJ*:-4 

3.6E-5 

l.lE-5 

3 .6E-6 

l.lE-6 

It 

5 . 0 E-' 

5 3E-5 

l,8E-5 

8.8E-6 

3 .4E-6 

l.lE-6 


-1.7C-3 

3.0E-4 

4.4E-5 

6.2E-6 

8.7E-7 

1.3E-7 

1 ^ 

1 lU- 5 

2.4E-4 

1.7E-4 

l,lE-4 

7.9E-5 

5.4E-5 

It 

G.3L-4 

1.9E-4 

3.9E-4 

7.4E-4 

1.7E-3 

3 .3E-3 


n Results for global IDC and IDDeC 






IDDeC 

iteration 

ntiniber 


p 

N^n 

IDC 

1 

2 

3 

4 

5 


4 

2.1E-1 

9.1E-3 

8.7E-4 

7.8E-4* 



1 

8 

8.0E-1 

1.5E-3 

2,4E-5 

8.9E-7 

3.8E-8 

7.6E-9 


12 

4.0E -2 

4,3E-4 

3,1E-6 

2.0E-7 

4 , 5E-9 

S.IE-IO 


16 

2.4E-2 

1.6E-4 

9 *88-6 

4.1E-6 

8.7E-7 

6.9E-7 


4 

4.1E-1 

3.9E-2 

4*3E"3 

8.9E-4 

6.3E-4 

6.0E-4 


8 

l.lE-1 

5.0E-3 

8.0E-4 

1.9E-4 

4.5E-5 

l.lE-5 

2 

12 

S.3E -2 

1.7E-3 

4.8E-4 

1 . 6E-4 

5 . 6E-5 

1.9E-5 


16 

3„0^,-2 

7.6E -4 

4.2E-4 

2.7E-4 

1.4E-4 

6.3E-5 



’ .n, 1 

2.7E- 2 

2 . 3E- 2 

5.3E-5 

2.7E-4 

3.0E-4 

) 

3 

^ .1i 1 

2.6E-3 

2 2E»4 

l.lE-5 

5.3E-7 

3 .OE-8 


1 2 

: .1 -2 

2.1E 3 

4.0E-4 

5.7E-5 

6*4E-6 

7.7E-7 


’ o 

2.2. -2 

2. JE-3 

4.1E-4 

1.4E-4 

6.8E-5 

6.6E-6 


f t 
i 

j .oP.-l 

7,8E-3 

1.6E-3 

3.6E-4 

7.7E-5 

1.6E-.5 

i 

) 

; 2 

n '"’71 n 

/ , 

5.4E-3 

3.9E-3 

2.8E-3 

2.0E-5 


\ f 

i 

2.0J-2 

7,0E-3 

1.4E-2 

2.9E-2 

5.CS-2 

9.7E-2 


g-j j'A; ' !'’ i'^ e sults for global IDC, and IDDeC 






IDDeC 

iteration 

number 


p 


IDC 

1 

2 

3 

4 

5 


4 

2.1E-1 

9.1E-3 

8.7E-4 

7.8E-4* 



1 

8 

8.0E-1 

1.5E-3 

2.4E-5 

8.9E-7 

3,8E-8 

7.6E-9 


12 

4.0E -2 

4.3E-4 

3aE-6 

2.0E-7 

4.5E-9 

8. IE- 10 


16 

2,4E-2 

1.6E-4 

9*8E-6 

4.1E-6 

8.7E^7 

6.9E-7 


4 

4.1E-1 

3.9E-2 

4*3E-3 

8.9E-4 

6.3E-4 

6.0E-4 


8 

l.lE-1 

5.0E-3 

8.0E-4 

1.9E-4 

4.5E-5 

l.lE-5 

2 

12 

8.3E-2 

1.7E-3 

4.8E-4 

1.6E-.4 

5.6E-5 

1.9E-5 


1 6 

3 . 0 1- 2 

/-6E -4 

4.2E-4 

2.7E-4 

1.4E-4 

6.3E-5 



‘ Ju 1 

2 

2.3E-2 

5.3E-5 

2.7E-4 

3 . OE-4 



•I , 3 1 

2 o 6£^ 3 

2 2E-4 

l.lE-5 

5.3E-7 

3.0E-8 


' 7 

.3 -2 

2. IE 3 

4.7E-4 

5.7E-5 

6.4E-5 

7.7E-7 



7.', -2 

2.3E-3 

4.1E-4 

1.-4E-4 

6.8E-5 

6.6E-6 


f * 

j * -1 

7„8E -3 

1.6E-3 

3.6E-4 

7.7E-5 

1. 6E-5 

. 

1 ^ 

7,1 : 2 

n " -p n 

/ m vJ 

5.4E-3 

3.9E-3 

2.8E-3 

2.0E-5 



-2 

7.013-3 

1.4E-2 

2.9E-2 

5.03-2 

9.7E-2 


[0,l] n,nd th^ values for U w^ro 4^ 8, 12 and 16. For tlis 
IDDoC, dogr-^'^ n vas t ken as N. Maximum a’bsolut'^ 

i^rrors ly n'othods were noted. Wp found that the IfCs 
With s -= 2 and 3 failed to improve thp results for th« IDG 
with s - 1 ( the raeaning of s as in the paragraph helow th”^ 
‘(luation (4.4.13)). 

h' fi'bl-'s 4. ^.10 and 4. 4.11 we display the errors 
nrofUu'i‘<l by th^ la^thods based on the local and global 
V rr.JO'io of IDG with s 1. 

*'ur Vv^sults show that deferred correction methods 
\ ■<.'"! ou th-’ locol approach arc pr'^f enable to those based 
on tb ’ ,'lob''l '•'pproach. They show also that the choice 
u -■ produc'jG worse results than p = 3 and that the choice 
')*' ! i '1' Vilu'^n such as 16 for m should be avoided. 

* ••’•3 r or linear Volterra Equati ons 

Ja this section we describe two deferred correction 
t -chni lu T ■’pplicabl^ to th^ nonlin-^ar Volterra integral 
'll a ’ tio 

X 

(..•..i;) y(x) - x,t,y( t) )dt = f( x), 0< x < a. 

0 

'I’th th ' us-^ of Gregory's formula (4.4.10) we obtain 


tl. ■ 'UT •'OXT'iate equation 



10 5 


r-1 

( 4 . 4 . 15 ) - h Zk(rh, 3 h,y.) - ^ hk(rh,0,Y ) 

X J u 

j=l 

- i hk(rh,rh,Y^) = 

where p is a chosen fixed constant and 

P 

(4.4.16) 0^ ^ = h S c^[ A®k(rh,0,Y^) + (-1)® ^ ®k( rh, rh, Y^ ) ] 

X jj^ p o U X 

s=l 

is the correction with the values of given in Section 4.1. 

There are two versions of deferred correction. One 
is the 'local' IDC which has been described by Kershaw [1974] , 
and the other version is the 'global' IDC which is introduced 
below. 

The local I DC 

¥e begin by finding the solution Y^°^ of equation 

( 4 . 4 . 15 ) with r = 1 and C =0. 

^ f P 

In the rth stage (2 < r < IT ) of the local IDC, we set 
the iteration count s to zero and perform the following four 
steps as many times as are necessary. 

Step I s With the use of values Y^°\ Y^°^, Y^^ \ 

Y^®^ and equation (4.4.16), compute 

X r y P 

Step 2: S-t Y^®^ = 



s tep 3 : 


C O » y 


With the aid of the values ^ Y^® ^ 

Y^®] and the new value, Y^®'^, obtain from 

( 4. 4.15) 


3te_p_4; Set s = b -!- 1 

For a giver, r, bhe sefiuencG of iterates Y^'^^,.,= . 

obtoiined by the above process is not guaranteed to converge, 
and so Y^,^^ is accepted as the final approximation to 
y( x^) . Iho accepted value is then denoted by Y^®\ 


The globa l IDG 

In the global IDC„ correction is done only after the 
soliation is found at fall the grid points by the basic 
discretisation nethodo The correction steps for this IDO 
'ire the same as those of the local IDO except that now each 
of ytopc 1 to 5 is perforiaed for all r = 2(l)IT. For i = 1 , 
2 and 3 lot us therefore define S tep i ' as "For r ::r 2 ( 1 . )K 
Step i of the local IDG ” , 

In the global IDC we begin by finding the solution 

/ \ 

r = 0(1 )H of ecLuation ( 4 . 4 . 15 ) with C = 0 ; this 
1 r y p 

is the sarao as using the trapezoid method for ( 4 . 4 , 14 ). 

Then, for s = 0, 1 , 2, ... do the three steps 
Step 1 ' to 3'. 

Since the seluonce of vectors y''®^, y(T)^ 
produced by the above IDC may not converge, we take 
as the final approximation, .and rename it as 



l.'ur rcric a l Results 

'ic; consider th(3 numerical solution of six problems 
by the local and global approaches of IDG and the corres- 
ponding IDleC mothodSo The first five problems are 
Problems 1 - 5 of Section 2,3-. 2 and the sixth problem is 
th.!; one given in Section 2,4. The range of integration 
’.r:is taken as [0,1], and the maximum absolute errors are 
tabula tod in Tables 4»4«12 to 4.4«19. 

In the tables >ro hove encountered so far j.n this 
thoeis, the degree m of polynomials used in the defect 
correction procedure has only been even and the results 
obtaino-d have been highly accurate. However we m.ay allow 
n to ta,ke odd numbers also in or-dcr to obtain equally good 
results. Table 4. 4,; 12 iD.lus-fcratos the good perfonaance of 
a glo’oal IDG bo-sed IlleC with odd values for m . 

In Table 4.4,, 13 we have listed the results of the 
globa?. irC with p — 1 and the corresponding IDDeG with 
m ~ 6 , 

Hg contrast the pei-fomancos of the global and local 
lie approaches and compare those of the corresponding- 
IDDoC nethods in Tables 4.4.3.4a and 4.4.14b. The local IDG 
yields more accurate results than the globa,! IDG. But the 
Ibl'oG baced on the global IDG produces results almost as 
good as the IDDcG baaed, on tho local IDG. It is interesting 
to observe that each successive iterations of IDDeC improve 

the previously obtained revsults in spite of the irilure of 
the iterotec of the underlying global IDG to converge. 



Table 4.4*12 Global IDC with p = 2 and IDDeC with odd degree polynomials 
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Table 4.4.14a Local IDC with p = 4 and IDDeC with 










•able 4,4.14b Global IDC with p = 4 and .iLbeC v/it:h 
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Table 4.4.15 d.ispla^’r:? the resiilts of IDDeC with m = 8 
o,nd the underlying local IDC method with a single deferred 
correction invol-ving upto fourth order differences in 
G-i* ? ■ s f 0 rmula , 

Table 4.4.16 shows the power of the IDDeC method in 
’producing highly accurate solutions for relatively large 
stepsizos. 

Table 4.4.17 reveals the fact that for Problems 2 
and 3; increase in the degree of polynomials gives rise to 
an increase in accuracy of results by the IDDeC method. 

Tables 4.4.18 and 4.4.19 show that the performance 
of IDDeC based on IDG(p,p) for p = 1, 2, 3 and 4 generally 
becomes hotter as p increases. 

For pra.ctical purposes for the numerical solution of 
Tolterra equations, ire prefer IDDeC method based on the 
local IDC to that based on the global IDC. Moreover, we 
choose tho parameters m, p and s of the IDDeC method such, 
that n = S or 3, p = 3 or 4 and s = 1 or 2 03aly. 



'"'artZ.o 

\.lb 

IDDeC v;iti: 

1 m = 8 and 

the underlying local 




IDC (4,1) 

method 







1 

IDC 

... . . 1 

IDDeC 

iteration 


i lOijlerri 

K 

(4,0) 

( 4|> 1 ) 

1 1 

[ 

3 

5 

1 

b 

8 . 2t';— 3 

j 

iK— 

4E-7 

3EZ-11 

7E-13 


1 u 

2.0H;_3 

2E-~5 

2E-8 

7E-13 

2E-13 
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I*"*!! 

2. ;..!E-3 

2E-5 

! 6E-8 

1 

9E-8* 
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7,08-4 

3E-6 

i 2E-10 
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1.7E-3 

IS- 4 

1 lE-7 

1 

IE- 11 

4BZ-12 
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4„3Fi-4 
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; 8E-10 

5E-14 

lE-14 


^■5 

2. IE- 3 

lE-5 

i 8E-'’ 

9E-10 

4S-12 
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5,3E,-4 

2S-6 
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! 2E-8 
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8E-14 
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3„6E-3 
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I''’"' ' ■ 
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1 2E-6 
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8E-12 
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8.9E-4 

3.S-7 

1 5E-8 

4E-12 

5E'- 14 
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Table 4.4.16 Performance of IDDeC for small values of 
m = N based on local IDC (2^2) 


iteration number 

Problem 

n 

1 

3 

5 

1 

4 

6 

2.3E-4 

4.9E-5 

3.5E-8 

1. 5E-8 

9,7E-12 

4.0E-12 

2 

4 

6 

1 . 6E-4* 

4.0E-S* 



3 

4 

6 

2.9E-6 

5.3E-8 

1.8E-6* 

2.9E-9-* 


4 

4 

6 

3 .6E-6 

8.2E~7 

1.6E-7^ 

7.2E-10* 



''-x 

2.4E-6 

9.4S-10 

3.2E-13 
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S 


5.5S-7 


1.9E-10 


5.9E-14 








i' abi e 4.4.13 Resxjlts for problem 1 by IDDeC with m = 6 


based on local IDC with p = i, 2>. 3 and 4 


1 "’^ 

n 
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IDDeC 


\ IP / J 
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(P,3) 
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5E»5 

6E-8 
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48 
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Table 4.4,19 

Results for problem 

6 by IDDeC 

with 

m = 6 

based on local IDC 

■with p = 1, 

2, 3, 

4 and 

K =12 




•id^=ebl'ern: 

p 

IDC 

(p.p) 

IDDeC 
(p. 1) 


IDDeC 
(p, 3) 

1 

8. OE-3 

1.4E-4 


2.1E-7 

2 

3 . 7E— 3 

8 . IE- 5 


2.1E-8 

3 

4.22-3 

4, 2E-6 


1 . 4E-S 

4 

4. 2S-3 

1.9E-6 


1 .3E-S 


4 . 5 C onclusion 

We have pnesenijed "the meihod of IlDeC "which, is 
another new iterative technique for improving the accuracy 
of nunerical solution of operator equations; we used the 
Fredholm integral equations to 8xp3.ain the ideas about the 
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method and analyze its asymptotic hehaviour. 

The IDDeC method is a generalization of the IDO and 
IDeC methods. It tias a faster convergence than its consti- 
tuent methods, hut is free from, their disadvantages. For 
instance, IDTeC removes the uncertainty in convergence of 
the iterates of IDC, and it is less expensive than IDeC in 
that the number of applications of the defect correction 
to obtain a given order of accuracy by the IDTeC is- less 
than that required by the IDeC on which the former method 
is based, 

¥e have applied IDDeC based on the IDC methodj that 
use quadrature rules for FIBs and VIEs, In the case of 
linear FIBs, we have considered the underlying IDC methods 
that are based on Simpson’s as well as the trapezoid rule 
with correction. But in the case of linear and nonlinear 
VISs, we have described two approaches of IDC based on 
G-regory's formula. 

Now we shall remark on the computational experiments 
on IDDeC methods. When applied to FIEs, IDDeC based on 
Cregory's fomula produced confirmatory results for our 
theoretical analysis on the order of convergence. It was 
also found that the above method is easier to implement 
and is more desirable than the IDDeC based on Simpson’s 
rule with correction. As for the application to linear 
VIEs, the IDDeC m.ethod based on the local IDC performs 



119 


better than the one based on the global IDG Howevery 
in the case of nonlinear VISs, the above two IDDeC methods 
Tji’Oduce dually good results in spite of the superiority 
of the local IDO over the global IDC , 



CHA]?DBR 5 


THE METHOD OP SUCGBSSIVE EXTRAPOLATED 
ITBRATBI) DEEBRHED DBEBCT GORRBCTIOH 


5.1 Introduction 

In Chapter 4 •we introduced the method of Iterated 
Deferred Defect Correction (IDDeC) -which is the result of 
applying the method of Iterated Defect Correction (iDeC) 
on the technique of Iterated Deferred Correction (iDC). 

¥e also found that when applied to Eredholm integral 
equations of the second kind, the IDDeC method based on the 
IDC that uses Gregory’s formula possesses an asymptotic 
expansion in even powers of the stepsize h. 

In the next section we introduce a new method which 
is a judicious combination of the three standard techniques 
for iterative improvement, namely deferred correction, 
defect correction and Richardson extrapolation. 

5 . 2 Su ccessive Bxtranolated IDDeC (SEIDDeC) 

We apply Richardson extrapolation on the method of 
IDDeC, and name the resulting technique as the method of 
Successive Extrapolated Iterated Deferred Defect Correction 
(SEIDDeC). 
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We now present this new method SBIDDeC as applied 
to Fredholm equations of the second kind. Start with a 
hasic stepsize of h = (h-a)/!)!. Find the nmerical solution 
of (2.2.1) by the ID0(p,p) method that uses the trapezoid 
rule with Gregory correction consisting of differences of 
order not greater than p and involves p iterations of 
deferred correction. Then apply the method of IDeC, with 
a fixed degree m of polynomials, on the above method of IDG, 
Denote by Y(x,h) the solution thus obtained. 

The order q of accuracy in this solution is given by 

(5.2.1) q = min[( s+l)(p+2), m+2], 

where s is the number of times defect correction has been 
performed (see Theorem 4.2 in Section 4.3.2), For a suit- 
able choice of s, the solution Y(x,h) is of order in4-2 and 
an asymptotic expansion for the global error Y(x,h) - y(x) 
is given ty 

(5.2.2) Y(x,h) - y(x) = C^^2^x)h^'^2 + x)h“'^^ + ... 

How apply Richardson extrapolation as in Section 3.3. 
Continue the extrapolation process until convergence, to 
within a user-prescribed tolerance, of the extrapolated 
values is obtained. If this is not achieved in a few 
iterations, say 2, the entire procedure is repeated with 
a smaller basic stepsize. 
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The asymptotic expansion for the method SSIDDeG is 
given hy 

(5,2,5) Y^(x,h) - y(x) = ^^^^+2+2 

+ •••••. 

•where j is the order of extrapolation. This is analogous 
to e^iuation (3.5.3) vhich gives the asymptotic expansion for 
the SEIDeC method. 

The SBIDDeC method described above is denoted by 
S'0IDDeCj^(p,s) . Ve observe that SSIDDeC is to IDDeC what 
SEIDeC is to IDeO . 

5 , 3 Apnlication to Predholm Bcuations 

In practical situa';ions we start with an IDC for 
which p is 2 or 3> and the'n apply IDeC for which m is 6 or 8, 
and the nimaber of iterations s is 1. In fact we can 
choose 2(p+2)> m4-2 or 

( 5.3.1) P> 2/2 - 1. 

This ensures that- the IDDeC(p,l, will be of order m+2 (see 
equation 4. 3 . 19a) ) . Purthermoro, Richardson extrapolation 
may then be applied on this IDDeC. The method described, in 
thp above paragraph is SBIDDeC^(p,l) for which p and m, 
satisfy the relation (5.3.1,. Since the defect correction 
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procedurs is costlier than deferred correction, we always 
try to minimize the number s of iterations to one. 

For our numerical experiments, we considered 
IDDeC(p,l) methods with m = 4 and 6 based on IDC(p,p) with 
p = 2 and 3 respectively and the corresponding SBIDI)eC^( p, l) 
methods. The results for IDDeC and SBIDDeC correspond to 
those of IDeO and SEIDDeC (see Tables 3.4.1 and 3.4.2) and 
they are identical except for the entries which have 
reached the machine accuracy. Hence we do not display here 
the computational results for IDDeC and SEIDDeC. 

5.4 Conclusion 

In this chapter we have presented the SEIDDeC 
technique which is a g^^neralization of the methods of 
deferred correction, defect correction and Richardson 
extrapolation. The order of convergence of this method is 
m+2+2j, where m is the degree of polynomials used in the 
defect correction procedure, and 3 is the order of extra- 
polation. 

in advantage of the method of SEIDDeC over the 
method of SEIDeC is that the former usually requires at 
most one iteration of defect corraction whereas the latter 
requires m/n iterations, n being the order of the basic 
discretization scheme. 
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Ths S,3IDDeG method is applicable to any basic dis- 
cretization scheme for the num,erical solution of an operator 
eluation vrhenever the scheme admits an asymptotic expansion 
for the IDDeC method on which the SSIDDeC is based. 


GHAPTPIR 6 


ThJ :i?;P:IO D OF ITBRATB I) DEFECT C ORRECTION OK aXTRAPOLAPlO I 


6,]. Intro f luo tio n 

In Chapter 3 we applied Richardson extrapolation on 
the method of IDeC and thereby introduced the method of 
Gucccssive 'Extrapolated IDeC (SBIDeG)^ In this chapter we 
combine the methods of extrapolation and IDeC in the reverse 
order, by performing IDeC on the method of Richardson extra- 
polation, The resulting technique is named as the method 
of Iterated Defect Correction on Extrapolation, This new 
method is presented as applied to the Fredholm integral 
equation 

b 

(6.1,1) y(x) -j’k(x,t) y(t) at =f(x), a< x< b. 

a 

For the above FIE, quadrature methods based on let^ton- 
Ootes quadrature formulae were discussed in Section 2,2, 
and the principle of Richardson extrapolation was explained 
in Section 3.3. The method of IDeC on Extrapolation is 
introduced in the next section, and its asymptotic expansion 
is derived in Section 6 . 3 . 
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h . 2 It-j ratod lef act C orre ction on Extra,polat i on (iDeOB) 


Considering extrapolation of the first order, we now 
I. r^sont th'-; method of Iterated Defect Correction on Extra- 
polation (IDcCE) as applied to a quadrature method Quad(h) 
v/ith a st-.'psize h for the numerical solution of (6,l„l), 


We first solve (6.1.1) by the methods Quad(h) and 

i"; \ ( 0 ) 

^.uadi-^i, and then obtain extrapolated values at the 

points x^ = a + ih, i = 0(l)N by Richardson extrapolation 
of ord'r one. 

As in Section 2.2,1, we proce-^d to apply the method 
of IDeC by constructing Ji polynomial function 
which intverpolates the points ( x^^, Zy^), i = 0(l)S. ¥e 
then solve the NP 


b 

(6,2.1) y(x) -j' k(x,t) y(t) dt 

a 

b 

= -/k(x,t)p^°>(t) at 

a 


first by duad(h) and then by '^uad(-|) , and finally obtain 
the resulting extrapolated values ^ at x^^, i = 0(l)K. 

Wo then compute an improved solution with the 

use of ?aadunaisky ' s principle 


( 6 . 2 . 2 ) 


,( 1 ) 7 ^°) 




The procedure above is a simple defect correction on 
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"• Ttrapolation. 

Ropeatad applicatioii of defact correction leads to 
the ethod of Iterated Defect Correction on Extrapolation 
( IDcCB), 


Practic al difficulties 

There are two difficulties that are encountered in 
the i.rplementation of the IDeCB method discussed above. 

One difficulty is caused by the nonavailability of extra- 
polated values at the midpoints of the original gridpoints. 
This difficulty arises when finding the solution of IIP by 
^uad(-^) which calls for a knowledge of the values of the 
function on the right hand side of (6.2.1) a.t all the 
in. 

= a + 1-2 , i = 0(1)2]!. Besides, we wish the values 

P^^Hx-') at the above points xj to have the same order of 
m i 1 

accuracy as the extrapolated values available at the original 
grid points. To avoid this difficulty, we first find 

'‘i^d2\o) at X.' , 3 = l(2)2N-l, using the equation 
2 3 

N 

(6.2.3) = £ w,k(x.,t )Z),°^ + f(x.), 

3 -1-3 11 3 

i=0 


which is obtained from (6.1.1)? we then construct T^°^(x) 
interpolating the 2N+1 points ( Xj; , Z^^°^), i = 0(l)]! and 



3 ^ 1(2)2^T-1. 
D 


1 


Another difficulty is that the IDeCE requires extra 
storage. Por an efficient computation of the solutions 
to Op r.-nd HP hy a quadrature method, the LU decomposition 
oP ft' Matrix of coefficients of the system of linear equa- 
tions rciBulting from the Op has to be saved, luring any 
dofect correction proceSvS, IDeCS thus requires a separate 
storage for a ( ZE+l) x (21^+1) coefficient matrix in 
connection with the use of Quad(-^^) in addition to the 
storo.g'-' needed for an (N+l) z (N+l) matrix for the use of 
Mun.d(h), This is not so in the case of the SSIDeC method. 

^ * 5 A gym-ptotic Error Sx-pansion 

Wo now give a derivation for an asymptotic expansion 
for the global error in the solution to (6.1.1) obtained 
bj’- the IDeCE method described in the last section. Since 
this derivation is similar to tho,t for the asymptotic 
expansion for the IDeC method discussed in Section 3-2, 
the details will be omitted. 


v/e assume that rn is an even integer greater than two. 
Let the errors in the solutions of CP (6.1,1) by the metho-ds 
4 uadfn(h) and '4uadm(v) based on the orapezoid rule, and the 


method of extrapolation be denoted by h and n; 


2 


,.(1) 

i, h 


respectively, and those of KP (6.2.1) by h and 

( 1 ) 


2 


e; . Then for the points x. 
i,h ^ 1 


a - ih we have 
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'^i,h “ A^(Xj^)h^ + ..., 


h = A 2 (x^).^ + A^(Xj;)^g + .,., and 


( 6.5.1) 




Consequently, 

(6.5.2) - y(x.) = + B^.^\x.)fa^ +' ... 

i 1 4 X n X 


Now th.3 relation between the errors by the Quad^(h) 
method applied to Op and NP is given by 


(6.5.3) 
Therefore, 

(6.5.4) 


e. , = -r. , + 0(h^ 

i,h i,h 


+ ••• 


and analogous to (6,5.2) have 

(6.5.5) = B^^^\x..)h^ + p^g^\x.)h^ + ... 


The use of Zadunaisky’s principle (6.2.2) yields 


(6.3.6) Z\-^^ - y(x^ 


^ tr ^ B^g^^(x^)h° + Bg^^(x,. )h® 4- ... 


Moreover, the required error expansion for the IDeCB 
method based on Quad^(h) is 

(6.5.7) - -y(x^) = B^®+^\x^)h'^ + B^®2^^(x.)h^'^^ +'.., 
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where q. = nin[2( s+l)+2, m+2], and s = 0, 1, 2, ... . 

In the case of the Qa.adg(h) method "based on Simpson’s 
rule, the asymptotic expansion for the corresponding IDeGS 
method is given by (6,3.7) with q. = niin[ 4( s+l)+2, ni+2]. 

• 4 C oncl u sion 

In this chapter we have presented the method of IDeO 
on Extrapolation as applied to the q.uadrature methods for 
the numerical solution of Fredholia integral eq, nations of 
the second kind. We have also derived asymptotic expansion 
for the above IleCS method. 

A common feature of the IDeCE and IDeC methods is 
that each defect correction increases the order of accuracy 
in their solutions by the order of the basic discretization 
scheme, and that the attaina"ble order is limited by the 
degree of polynomials, 

When the methods SBIDeC and IDeOE are compared, 

SEIDeC is found superior to IleCB, The disadvantages of 
IDeCS are that it requires more work and more storage than 
SBIbeO, and that its order is limited. 

To summarize, we find that the order in which the 
methods IDeO and Richardson extrapolation are combined 
does matter and that the order which is used in SBIDeC is 
to be preferred, . 


CHAPTER 7 


THE MBTHOP OP SUCCESS lYE UPDATED ITERATED 
EBFEGT COREBCTIOR 


7.1 Introduction 

In this chapter we present a new techniq,ue called the 
method of Successi-ve Updated Iterated Defect Correction 
(SUIDeC). This technique is a variant of the relatively 
new IDeC method. The SUIDeC method makes use of inter- 
polatory cubic splines instead of polynomials which have 
been employed in the defebt correction procedures discussed 
in the previous chapters. 

V/e introduce the new method as applied to the follo’"- 
ing nonlinear Yolterra integral equation of the second 
kind: 

X 

(7.1.1) y(x) = f(x) +j' k(x,t,y(t)) dt, 0< x < a. 

0 

As in Section 2.3, we denote the unique solution of this 
equation by y(x). 

Let N be a positive even integer, h = a/U and 

= rh, r = 0(1)N. Let us represent an approxima- 
tion to y(Xj,) by and denote f(x^) by f^ and k(x^,t^,T^) 
by k ■. Then loble’s method is given by the equations 

T y 3 
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X 

( 7.2.1) y(x) = f(x) + j X, t,y( t) ) dt, ^ 

0 

¥e then compute a solution value at for 

this e<luation by Foble’s method, and its iterative improve- 
ments Y^^^Cs > 2) by the method of IDeC. ¥e denote by Y^ 
the improved value obtained at the end of the final itera- 
tion of the IDeG. 


■S ka i-th . stage 

■Having determined the values Y^ at for j = 0(1) i-1, 
we use (2.3.4a) or (2.3.4b) of Noble’s method as required, 
and find Y^^^ as an approximation to y(x^). To improve 
tho solution value just found, we estimate its unknown 
error y( x. ) - yJ. V/e fit the points (x., Y.), 

1 J. J J 

1 = 0(l)i-l and ( x . , Y^ by an interpolatory cubic spline 
function x) , and construct the problem 


(7.2.2) y(x) = - /k(x, t,p\^'^^( t) ) dt 


( 1 ), 


X 


+ j‘k( X, t,y( t) ) di 
o 


0 < X <_ x^. 


X 

Since P^;^\x) - rk( x, t,P t) ) dt - f ( x) is small, 

0 

(7,2.2) is a neighbouring problem (NP ) of (7.2.1). The 
exact solution of this IP is ohe known function p\ -ix). 
Now, starting with updated values Y^, j = 0(l)i— 1, and 



using ( 2 , 3* 43*) or (2, 3# 4b) of Hoble^s msthod as rsquired, 
we find th- solution of Fp (7.2.2) at llie error 

in this solution value is known, and is equal to 

J-i vXj.; - r\^ , 


Since both (7.2,1) and (7,2.2) are solved by the same 
ihthod, we postulate that the error y(x^) - is 

approximately equal to the known error P^^\x.) - 

U X X 

Therefore, 


y^2) ^ y(^l) ^ (p(.l)(x.) - r^l)) 


is a better approximation to y(x.) than y [^\ 


The above technique may be repeated to compute 

(s+lV 

improved- approximations' Y^ ' ' f or s - 2, 3, , . . , "by ' , 

performing the following 'steps; - ; ■* ^ • 

(i) Construct the FP 

■ ' - 

(7.2,3) y(x) =: x)- - /k(^x dt 

... ■ 

.. .» •+*/ (x, t, y( t) ) dt,. /O 

0 

'( q \ ' ' ' " ; 

where P- '^(x) is an interpolatory cubic spline function 
fitting the points (x,, Y^), j = 0(l)i-l, and ( x^, Y^^O. 
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{ ii) Find the solution value ^sing the 

relevant equation of Hoble's method. 

(iii) Compute the improved approximation 

(7.2.4) y^s+1) ^ yU) ^ 

For a preassigned tolerance TOL, -we denote by s^^ the 

value of s for which j ^ satisfied. 

¥0 then take updated solution value at 

for ( 7 . 2 . 1 ). 

Che values i = 0(l)K produced by the method of 
SUIDeC form the numerical solution of (7.1.1). 

Im-plementation reraarks 

In our implementation of SUIDeC, we calculate the 
end-derivates required for the cubic splines Hx) of 
(7.2.3) by making use of at least four ordinates (see 
equation (2.11) in Tewarson[l980 ] ). Hence we need three 
starting values Y^ and Y^ for SUIDeC. 

x 

The integral J k( X, t,P..®^( t) ) dt in (7.2.2) and (7.2.3) 
0 

is evaluated using standard routines such as those provided 
by NAG. 
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Post 

We shall now find the major difference in the cost 

of the global IleC method that uses splines and the SUIDeO 
method. 

In the case of the global approach of IDeC using 
Splines, an UP is of the form: 

X X 

(7.2.5) y(x) =P(x) -/ k( X, t,P( t) ) dt + J k( X, t,y( t ) ) dt, 

0 o 

0 < X < a 

where P(x) is a cubic spline function interpolating the 
points ( Xj^, YjL)> i = 0( 1)N and the in turn constitute 
the most accxirate solution available during the defect 
correction process. Thus during any iteration of the above 
IDeC only one spline function (through F+l points) is to be 
computed. On the other hand, in the case of SUIDeC an NP 
is constructed for each subinterval [0, Xj^], 4 < i< Y, 
and this NP requires a spline through i+1 points. Thus 
for each defect correction the number of splines to be 
computed in SUIDeC is about whereas only a single 
spline is to be computed in the global IDeC method. 

Apart from this difference, the costs of the methods 
IDeC and SUIDeC arc essentially the same. 
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"^•5 Ap-plication to Volterra Squations 

¥3 applied Noble’s method, the global approach of 
IDeC using splines, and our method SUIDeC on four problems; 
the first three are Problems 1 - 3 of Section 2 . 3.3 ■ : 

(Netravali[l973] and Baker[1977] ), and the remaining one 
is Problem 4 , from Baker [1977], given below; 

Problem 4 

X 

y(x) = 2 - + /e^‘”^y2(t) dt, ^1^11- 

0 I 

! 

exact solution; y(x) = 1 < 

Problems 1 to 4 were solved for N = 16 and 32 and the I 

I 

absolute values of the errors n = 3(1)N ' 

were computed. In Table 7.3.1 we give ft ese values ' 

corresponding to a few points in [0, l]. The value chosen 

for the tolerance TOl in the SUIPeC method was 10“^^. The I 

word 'iter' denotes the number of iterative improvements i 

done to improve the solution values. 

, |i 

7 . 4 Oonclusion ; 

’/o have described a version of the method of Successive I 
Updated Iterated Defect Correction as applied to (7.1.1).. ; 

We compared the results with those obtained by the method 
of IDeC involving the global approach and the use of cubic 
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splines, and found that the local approach (SUIDeC) yields 
bettor results than the global. 

Comparing the performances of the three methods, 
namely IBeC with polynomials discussed in Section 2.3, its 
analogous IBeC with splines, and SUIDeC of this chapter, we 
find that the first IBeC is the best and the cheapest. 



Table 7.3.1 i=ibsolute erirors in the solutions by the methods 
of Noble, global IDeC and SUIDeC . 




Problem 

N 

X 

Noble 

IDeC 

iter 

SUIDeC 

iter 



.5 

.l6(-8) 

.50(-ll) 


.ll(-iio) 

2 


16 

1.0 

.90(-6) 

.18(-6) 

3 

.23(-7) 

4 

1 

32 

.5 

.98(-10) 

.34(-12) 

3 

.75(-.12) 

2 



1.0 

.55(-7) 

.ll{-7) 


.16(-8) 

3 



.5 

.90(-7) 

.50(-8) 


.36(-8) 

3 


16 

1.0 

„ll(-5) 

.14(-6) 

2 

.57(-7) 

3 

2 

32 

.5 

.57(-8) 

.38(-9) 

3 

.30(-9) 

2 



1.0 

.66(-7) 

.97(-8) 


.49(-8) 

2 



.5 

.66(~7) 

.36(-8) 


.12(-8) 

3 


16 

1.0 

.26(-6) 

.73 (-7) 

3 

.77(-8) 

3 

3 

32 

.5 

.41 (-8) 

.2l(-9) 

3 

.64(-l0) 

2 



1.0 

.l6(-7) 

.44(~8) 


.45(-9) 

2 



,5 

.ll(-6) 

, .38(-7) 


.67(-8) 

3 


16 

1.0 

.63 (-6) 

,43 (-6) 

3 

.53 (-7) 

3 

4 

32 

.5 

.68 (-8) 

.23 (-8) 

3 

.12(--8) 

2 



1.0 

.38(-~7) 

.26(-7) 


.68 (-8) 

3 





GHAPTER 8 


RBLATIinB PERPOBMANQES OE METHODS OE ITERATIVE IMPROTSMENT 

Various iterative. Improvement techniques for soiying 
integral eqimtions are compared in this chapter, 

8,1 Applicability to Fredholm Integral Equations 

In this section we compare the following methods of 
iterative improvement which are applicable to Fredholm' 
integral equations: Extrapolation, IDG, IDeC, ItffieC, SBIDeC 

and SEIDDeC. 

The methods IDeCB and SUIDeC are not in the above 
list because IDeCE has already been compared with SBIDeC 
(see Chapter 6), and SUIDeC is not applicable to EIBs. 

We compare the methods with regard to the increase 
in order of convergence per improvement or correctionj the 
sequence of orders of convergence of the iterates, and the 
maxjjntJm possible order of convergence. We assuma that n 
is the order of the basic quadrature method on which 
techniques of iterative improvement are applied, and that, 
in the case of deferred correction methods, p is the order 
up to which differences are Included as correction to the 
basic quadrature formula used. ¥e summarize our findings 


in Table 8.1. 
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Table 8.1 CornTparlson of Iterati-ve Methods of Im-proTement 


method 

increase 
in order 

sequence of orders of 
iterates 

B xtrapola t ion 

2 

n+2, n+4, . 

IDC 

1 

n+1, n+2, p+2. 

IDeC 

n 

2n, 3n, . . . , 

IDDeC 

p+2 

2(p+2), 3(p+2), m+2. 

BBIDeC 

2 

m+4, m+6, ... . 

SB IDDeC 

2 

m+4, m+6, ... , 


As noted in the “able, if the iterative mithods for 
PIB(2.2.1) are based on the trapezoid rule, a siiigle 
Iteration of IDC increases the order of accuracy by one; 
that of IheC, by two; and that of II)DeC, .by p+2. 

I 

I 

Bach of Our methods IDDeC, SBIDeC and SBIDDeO is 

f: 

superior to Lindberg's technique, which improves the 
order of its iterates by n only. Moreover, the results 
obtained by our methods are far better than those by his. i 

' ' I 

(Contrast our results, for instance, in Table 4.4.4a with 
his in p. 79 of Lindberg.1976] . ) 

' is 

We recall that IDDeC is based on an IDG which is of 
order p+2, and that SB ID eO and SB IDDeC are based Respectively 
on an IDeC and IDDeC, each of which has the maximum order 
13+2. This order is attained after (m+2)/n-l applications | 
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of defect correction in the case of SSIDeC, and after 
(m+2)/(p+2) - 1 applications in the case of SBIBDeC. 

Thus, for a suitable choice of p, IBDeC needs only one or 
two applications of the defect correction process to 
produce results with error 

It is possible to modify our SBIBeC and SBIDDeC 
methods so that their underlying defect correction t echnique 
has a variable order of convergence, namely (s+i)n or 
(s+l)(p+2), instead of the constant order m+2. 

Bor general ptirpose routines for solving FlBs, it 
seems best to use IBDeC Or SBIDDeC based on the irapezoid 
rule. 

8. 2 Annlicability to Yolterra Integral Equations 

In this section, we compare accleration techniques 

applicable to YIEs,, 

Now that the SUIBeC method of the previous chapter 
has already been found inferior to the IDeC method of 
Chapter 2, we confine our attention only to the three 
methods IBC, IDeC and IDDeC. 

Using 9-11 significant digits, two problems in 
linear VIBs were solved by the three methods based on the 
trapezoid rule. Problem 1 is given below and Problem 2 is 
Problem 5 of Section 2.3.2. 
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X 

y(x) + y(t) dt « 2x , 0 < x < 1 

0 

exact solution; y(x) = 2x - x^. 

The maximum absolute errors obtained are given in Table 8.2. 
Ihe performance of IDDeC(2,l) is better than that of IDeC(2) 
hnd IDC. This is also ttue for the problem considered in 
iection 2,4 (see Tables ^,4.1 and 4>4.10), 

In the case of linear and nonlinear VIBs, Our 
Somputational results shOw that, to produce solutions of a 
given accuracy, IDEeC generally requires fewer applications 
Of defect correction than the corresponding IDeCi Thus it 

is cheaper to use IDDeC than IDeC. 

8.3 Final Remarks 

The method is IDDeC is a promising method with 
applicability to differential as well as integral equations. 
In the case of linear PIBs with smooth kernels and func- 
tions, SBIDDeC surpasses IDDeC. 

The hierarchical order of the methods as regards 
their capability to produce highly accurate solutions is 
as IDDeC, IDeC and IDC. 








T^le &m2 Maximum Errors by IDeC^ Local TDC, and IDEeC 






GHAPT3R 9 

APPLICA [i:iO]\l O F IDe O TO P ARTIAL DIFFPRUKT IAL JSOUATIOI^ 

9 f i Inuroduc -fc x on 

Vo no^’- turn our a b Lon bon co ap plica bn ons ol tlie 
H'-’thod of xt'^rabed d'^f^ct coi-X'^ction on tna Ixnite dxffc- 
■^once methods for partxr.l dxff e'^onl^al eo.uabxons, 

Zxduraxsky (1976) applxcd b .13 technique to partxal 
dxfferonbial oo^uatxjns foi’ scIp purpose of b^xia bing 
the glob'll G' rcr. Ix b' xs chap bor ue ^tbenpt to riodrxs^ 
bis method ond to us^ 7'' = reoulbi'ng algor 3 xhm xceratxv-lxp 
the results of thxr chapter o,’’? of an experimeabjl 
na Lurcc 

In tb^ next cecixcn ue pxpJa.'n una modxf xca .xon i" 
Zadunaxskyc prxncxpl? as applxod to p^irabolxo parxxai 
dxfferentxal oq.uatxons„ by x icTocucing certaxn funct .ons 
called 'x-splxnes' and ' b-spl : nes ’ , 

Secbxops 9»3 o.nd deal yit^h the appi xcapxon o. 

ID'^C r ^specxxvely to nyp^^ibolxc and elliptxc partxnl 
dxfferentxal eciuabions. In case of ellxpbxc equabiOnSj 

we consxd^r IDeC wxth xh'^ axd of pxf^cewxse xnterpol a lore 
polynornxols instf^ad cf th'= r-splxn^s arid the t-splxxes used 
xn ID'=G for the oth-^r two typ-^s of partxal dxf f er'^nbxal 
siuatronc , 
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Computational examples illustrate the performance 
of the modified IDeC. 

9 . 2 Po rabolic Partial Differential Equations 

¥e consider the following initial boundary value 
problem (lETP) of parabolic partial differential equations 

u^ = f( X, t,u,u^,u^^), a<x<b, 0<t<T 

(9.2.1) u(x,0) = gp(x), a < X < b 

u(a,t) = g 2 {t), 0 1 t < T 
u(b,t) = gj(t), 0 < t < T. 

let us denote the unique solution of (9.2.1) by 

u( X, t ) . 

In order to introduce the discretization of (9.2.1) 
by the method of backward difference (BIM), we choose 
positive integers and M, let h = (b - a)/K and k = 
and define a grid G as 

(9.2.2) G = {(x^.'t^):x^ = a + ih, t^ = jk; i = 0(l)B, 

3 = 0(1)M}. 

Purth'^rraore , we corap uto tho approximate solution 
values u^ at the grid-points (^^>’^ 2 + 1 ^ using the 
finite difference equations: 
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= 0, i = 1(1)N-1, j = 0(1) M-1, 


with the initial and houndary conditions 

( 9.2.4) 

^ ^>3 ~ ^ ~ 0( 1)M. 

Zadunaisky (1976) extimates the resulting global 
error u(x ,t ) - u .5 iii the following manner. A bivariate 
interpolating function P = P(x,t) is found such that 
P(x ,t ) = u . Then the NP corresponding to the OP 
( 9.2.1) IS taken as 

( 9 . 2 . 5 ) = f(x, t,u,u^,u^^) H- P.^ - f( X, t,P,P^,P^^) 

with the samf^ initial and boundary conditions for the OP. 

The numerical solution v . to th - NP (9.2.5) is then 

■^9 J 

foun<3 "by th.-^ same BIM, using the eq.uations: 
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^1. .1+1 ~ ^1. n 

k 


(9.2.6) 


_ ffx t u „. ^4- , ? -.»,„J:t ^ > 

^3+1’ 1,0 + 1’ 2h 


^+ 1 . 0+1 " 


+ u 


p 

1 = 1(1)K-1, 3 = 0(1)M-1. 


Since th^ exact solution to NP (9.2.5) is P(x,t), the 
global error -V. . •^ P(x , t ) in the numerical solution 

1, J 1 J 

V IS known. This known error is taken as the required 
i» 3 

ostinate for the global »rror. 

’/a observ'^ thet in the abowe in'^thod of Padunaisky 
an explicit computation of th'^ bivariate interpolating 
function P(x,t) can be avoided provided that approximate 
values for the necessary partial derivatives of P are 
determined only at thc' points of the grid G-. Accordingly 
we detormino univariate interpolating functions which are 
named as th» 'x-splines’ and 't-splines'. An x-spline is 
defined to be a spline function of x which interpolates the 
dnta points (x , u ), i = 0(1)11 for which j is a fixed 
integer. Si-ilarly, a t-spline is defined as a spline 
iunction of t which int^rpoletes the points (t ,u ), 

3 = 0(1)m for which i is a fix^^d int'^ger. In order to 
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de'teriinn'^ 3.pproxiiTi?to values for the derivatives of P 
appearing in (9.2,6), we nerely construct M x-splines and 
K— 3 tir-splinos, and take the corresponding derivatives of 
th-^so splines. 

Ve have introduced the above rodification in the 
iterated defect correction rie#iod for partial differential 
■equations. 

In the case of IBVP (9,2.1) if M is greater than 30, 
we first irprove the solutions in the sub-grid consisting 
of Vq( < 30 ) tir^e-steps by an application of IDeC n-^thod, 
and then proceed further in a sirilar way, raking use of 
the improved solution for the niQth time step as the 
starting values for the sub-grid consisting of the next 
tine-steps. After obtaining improved values in the second 
sub-grid by another application of IDeC, we proceed to 
improve th'= solution values in the next sub- grid till the 
solutj.on values for Mth time-step are improved. 

Numerical results 

In our Implementation of IDeC the x- splines and 
t-splines were taken to be cubic splines (see lGwarson[l980 ] ) , 
W d'^note by IDpC(i) the IDeC m=‘thod with ’ 1 ' applications 
of Zaduiaisky's principle. 

Th? first fiv^' of the following six test problems 
are from Lindberg (1976), and Probl'^m 6 from Greenspan( 1974) . 
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Proble m s 1 to 3 

U, = U 0<X<TC, t>0 

w X.X 

u(0, t)=0=u('ti;,t) t>0 


with th- following initial vilue function g^(x)' 


(^) 


g^(x) = sir)(x) 

-t 


0 < X < TC 


jxact solution: u(x,t) = e ^sin x. 


(b) 


g^(x) = tin - x) ’^1^5,'^ 


exact solution: u(x,t) = £ (2n - l) 3 (2n l) (Pn-l)x 


n=i 


X 


0 < X < Tt/2 


( c) 


exact solution; 


g]_(x) =-; 


^ TC - X %/2 < X < 


% 


=4 Z 1 -nn 

" n=l(2) ^ 


sxn nx 


n 


P roblem 4 

u^ = (2 + sin T^xx^^xx + 0 < x< -rt/2, 

t > 0 

u(0,t) =0, u(|,t) = g“'^ t > 0 

u(x,0) = sin X 0 < X ^%/2 

‘^xact solution: u(x,t) = e ^sin x. 
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P robl' r ' § 

Burger’s eq.uation 

= -u u^ + 0<x<TC, t>0 

u(0,t)=u(n:,t) = 0 t> 0 

u( X, 0 ) = sin X ^ 1 ^ 1 

exact solution; 

oo 2 

4 L e”^ Ci-)sin(nx) 
u( X, t) = — 

OO y. 2 4- 

Io(i) + 2x e~ (i) cos(nx) 

° n=l ^ 

where the are th-^ modified Bessel functions. 

Probleri 6 

u. = u - x(t-l)e"^ 0<x<l, t>0 

u X-X 

u(x,0)=0 

u(0,t) = 0, u(l,t) = t t > 0 

exact solution: u(x, t) = xte 

Tabl"' 9.2.1 gives th-^* naxiiiur^ absolute errors in 
th^ nu]"’ 'rical solutions obtained when Problens 1 to 3 were 
solved by the backward difference method (BDM) as well as 
the IBi,G I Lothod with k = .05 and h = tc/10 and 7 t/ 20 . The 
number of points taken for th'^ determination of a t-splin^^ 
was 20 only. 



Maxumiam Error? by BDM and IDeC with k = . 05 


h = 7t/io 


h = 'V20 


1 roblem 1 


t 

RDM 

IDeC(l) 

BDM 

IDeC(l) 

05 

] .5E-3 

3.0E-4 

1.2E-3 

3.5E-5 

25 

6.3E-3 

1.3E-3 

5.1E-3 

1.6E-4 

45 

9.2E-3 

1.9E-3 

7.5E-3 

2.5E-4 

65 

l.lE-2 

2.3E-3 

9 .OE -.3 

3 .2E-4 

85 

l,2E-2 

2.5E-3 

9 . 6E— 3 

3 .7E-4 


’roblem 2 


Problem 3 


.05 

] .2E-2 

1.2E-3 

9.9E-3 

2.8E-3 

„25 

1.5E-2 

3.0E-3 

1.3E-2 

8.6E-4 

.45 

2.1E-2 

4.2E-3 

1.7E-2 

5.6E-4 

,65 

2,7E-2 

5.4E-3 

2.2E-2 

7.1E-4 

.85 

3.0E-2 

6,3E-3 

2.4E-2 

8.9E-4 

.05 

6.9E-2 

2,0E-2 

4. IE- 2 

1.2E-2 

.25 

3,0E-2 

7.3E-3 

1.8E-2 

2.3E-3 

.45 

2.3E-2 

4.2E-3 

1.4E-2 

1.2E-3 

.65 

2. IE- 2 

2.3E-3 

1.4E-2 

8.4E-4 

.85 

2.0E-2 

l.lE-3 

1.4E-2 

5.6E-4 
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The i"prov3i f'nt in the solutions by the IDeO is 

vj_rite cood for each of thi three prohlers solved. It is 

rot worthy that Zadunaisky's principl'= of prror estimation 

’'■'■■-n ■'/orkod well for Problem 3 for which thp initial value 

'‘’notion IS not smooth whereas th) method of Lindberg 

(1976) I'ail'c to estimate the errors reliably. We also 
i . o tiie results of IPeC(l) wiin n - a/lO are 
1 'utpi' than those of BDM with h = -n/ZO. However, IDpC(i) 

f or 1 ^ 2 did not yi'^ld better r^-sults than IDeC(l). 

The backward difference method was employed also 
to solv Problems 4 and 5 . For each time-step, the BPM 
7i''ldod a sy.eten of iion-linear equations which was solved 
by r-'TTton’s rtothod. Th'^ initial approximation for th^ 

;i' Hiod T,ras takc/i as th^ solution at th^ previous tirn-step, 
amd th iterations were terminated when the maximum change 
in the values of th® solution was less than a tolerance, 
/TOL . 

Tor each time-step, the system of non-linear 
equations arising out of the finite difference equations 
for tlip n ’’imhbour ing problem was again solv'^'d by Fewton’s 
m'^thod. The initial guess was taken as tht nu"’erical 
solution of the original probleri and iterations wer^ 
tcn-'ine b-’cl as in th'= case of tb^ original problem. 

Against each of th^ sal'^ct-’d timp-st'^ps in T^bl-^s 
9.2.2 and 9.2.3 are giv'-n the maxi’^'Un erroi in 

th-' solution of Problems 4 and 5 respectively and the 
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nrror t th^ central grid-pOint ; also are given the nunber 
cf iti-'^rations required by Fewton's method to solve the 
origin'’! problem and th'^ nupibr>r required -when solving a 
i''’ighbourin''' probl'^’a. Thf’ total numb'^'r m of th*^ tine-steps 
in a sub-grid for an application of IDeC was chosen as 10 « 

At the nid-points of the tne-lines, IDeC(2) produces 
brtter r^'^' suits than IDeC(l). Howev'^r, when accuracy over 
the •’ntiro grid is desirf='d th'^ IDeC(l) method is preferable. 

Th'' solution valu'-^s at points on x = -^ for Problem 6 
by the methods BDM and II)eC(3) with h = .05 and k = .1 and 
.01 arc given in Table 9.2.4. 

Ve notice that th'^ IDeC(5) m'^'thod with k = ,1 gives 
nuch bettor results than the BDM with k = ,01. As for the 
cost, tho number of arithmetic operations required by the 
IDeC(5) with k = .1 was a little less than the nmber by 
th® BDM with k = .01. 

We also solved Problem 6 by the BDM and IDeO methods 
for h = .0 5 and a few values for k. Th'^ errors in the 
computed solutions were measured by the Euclid ®,n norm in 
the grid comprising the points (x,t) with x = 0 (. 05 )l and 
t = 0( .1)1. The results in Table 9.2.5 indicate that in 
this case our iriiplementation of the IDeO method is quite 


eff 1 C lent . 



Tdble 9.2.2 Results for Problem 4 with h = 'K/20, k = .1 
and XTOL = .0001 



errors over time-line 



errors at 

centre 

hime-step 

BDM 

No. of 

IDeC(l) Iterations 
OP NP 

BDM 

IDeC(l) 

1 

3.2E-3 

9.0E-4 

7 

4 

2.9E-3 

2.9E-4 

5 

6.5E-3 

2.3E-3 

8 

6 

6.3E-3 

4.9E-4 

10 

4.2E-3 

5.8E-4 

6 

4 

4,2E-3 

5.5E-4 

15 

2.3E-3 

6.5E-4 

8 

5 

2.3E-3 

3 . 6E-4 

20 

1.3E-3 

1.9E-4 

8 

4 

1.3E-3 

9.5E-5 

25 

6.6E-4 

2.1E-4 

7 

4 

6.6E-4 

1.2E-4 

30 

4.4E-4 

9.0E-5 

6 

3 

4.4E-4 

8. 6E-6 

35 

2.8E-4 

9.6E-5 

6 

3 

2.8E-4 

7.7E-5 


\ib 1 e 'J . 2 . 3 Pesults for Problem 5 with h = k = .1 

and XTOL = .0001 


errors over 

time- lane 


errors 

at centre 

t i me- step BDM 

IDeC(l) 

No. of 
Iterations 

BDM 

IDeC(i: 



OP NP 




1 

1 . IE- 2 

3 .lE-3 

5 

2.0E-? 

3 .9E-3 

in 

2.2E-2 

2.1E-3 

35 

2.0E-2 

1.4E-3 

20 

1.6E-2 

9 . 2E-4 

25 

1.2E-2 

6.8E-4 

30 

8.9E-3 

4-9E-4 

40 

4.5E-3 

2.5E-4 

55 

1 .4E-3 

8. lE-5 


3 

CM 

4.9E-3 

6-5E-4 

3 

2 

1.9E-2 

3.3E-5 

2 

2 

2.2E-2 

6.4E-4 

to 

2 

2.0E-2 

7.9E-4 

2 

2 

1.6E-2 

7.8E-4 

CM 

2 

1. 2E-2 

6.4E-4 

2 

2 

8.9E-3 

4.8E-4 

2 

2 

4.5E-3 

2-5E-4 

2 

1 

1.4E-3 

8.1E-5 



Icble 9,2.4 Problem 6. bolu_ion values when x = 'i . 


t 

BDM 

h = .1 

BDM 

k = .01 

IDeC (3) 

k = .1 

IDeC (3) 

k = .01 

exact 

solution 

. ? 

.07808528 

.08143904 

.08178577 

.08187299 

.08187307 

.4 

.13029526 

.13367796 

.13405136 

.13406401 

.13406401 

. 6 

.16166281 

.16435400 

.16464544 

.16464349 

.16464349 

.3 

.17754286 

.17952423 

.17973062 

.17973159 

.17973159 

1.0 

.18239586 

,18379544 

. 18391273 

. 18393972 

.18393972 


T able 9.2.5 Errors obtained by the 5DM and IDeC methods withh=.05 


k 

BDM 

IDeC(l) 

IDeC (2) 

IDeC (3: 

0.1000 

0.6196 

0.1648 

0.0867 

0.0678 

0,0500 

0.2911 

0.0380 

0.0073 

0.0047 

0,0100 

0.1993 

0.0202 

0.0045 

0.0024 

0.0025 

0.1712 

0-0191 

0.0033 

0.0031 

0,0010 

0.0631 

NC 

NC 

NC 
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9 • L y,^~ ' rbolic Partial Diffarontial l i quations 

First, we consider la''- initinl boundary value problev 

u. =cu c>0, 0<y<l, 0<t<T 

u(x,0) = g^Cx) ® 1 X 1 1 

u(l,t) = g2(s:) 0 < t < T. 

In the notation of S-'ction 9.2, th-^ ’loap-frog’ 
ethod for (9.3.1) IS given by th ) -equations 

^'^i.n+1 ~ ^?-9.2 v1 _ Q ^±liX, Irl t-d = 0 

2k 2h 

subject to the conditions 

“x,o=sfa)> i = 0(l)E 

.=0(1)K. 

be consider th^ following conputationol exanple of 
non-linoar problem 

Probl'^m 1 

0 <:. < 1, t > 0 
u(x, o)=l-x 0 <X<1 

u(l,t)=0 t>0 

■'xact solution: u(x,t) = (l — x)/(l + t). 



w ' solv:d Problem 1 by the leap-frog raetbod 
IbeC method with h = .05 and k = . 015 . The bound* 

Vi /•» 1**^ 

0.1 T = 0 w'^r'^ taken as u(0,t) = 1/(1 + t), and 

wore calculated fro^i th'^ Lax— "'■'fendroff scho 


“l,l = “l,0 + 


+ « '^1+1,0 - 2“x,0 


The naxiiiun absolute errors for the t'^nth "t^. 
b.y the leap-frog, IDeC(l) and IDciC(2) methods ware 
9. 1^-6 and 4.73-6 respectively. 

Next, we consider the IBVP; 


( 9.3.2) 


0 < X < a, t > 0 
u(x,0) = f;]_(x) 0 < X < a 

u^(x,0) = f 2 (x) 0 < X < a 

u(0, t) = g 2 ^( t) t > 0 

u( a, t) = § 2 ^ ^ ^ • 


u — " u. I 

zx tt 


With the notation of Section 9.2, th'= nunerir. 

»n u. to the above IB' 

1, D 

thre-~-level implicit’ scheme: 


solution u. to the above IB'VP may be found by th 
3 


1 


h 




2'- '''°i “i, 3+1 


+ ( 1 - 


2 k) of u. + 

1 1, 3 


^ t 2 


^ 2 u , 
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wh'^ro A IS a 'r slaxation' factor and 
operator 


h'f VL .IS the 

i 1, D 


5 


. U. . e U , , - 2u . 4- U, 

1 1, 0 1+1, 0 i, 3 i-1, 3 


= the above implicit sch<^me becomes 


For our choice of ^ 
ntoble for all choices of h and k and requires the solution 
of a tridiaponal system of equations, for each tire-step. 
(For details, see AmGs[19773 and GreenSpan[1974] ) . 


Problems 2 to 4 w^re solved by the implicit scheme 
nentioned above withihe necessary modifications, as in 
pp 174-175 of Greenspan' (1974), for Problem 4. 


Problem 2 

^xx “ ^tt 

u(x, 0) = 2x + 0^, u^(x,0) = -e^ 

u(0,t) = e”"^, u(l,t) = 2 + ^ 

exact solutionJ u(x, t) = 2x + e . 

Problem 3 

^xx = ^tt 

u(x,0) = sin itx , u^(x,0) = 0 

u(0,t) =0 = u(l,t) 


°xact solutiont 


u(x,t) = sin iix oos lit. 



161 



In T'=i'bl? 9.3.1 W 2 give the errors, by the Hiclidean 
nort'i, in th’ computed solutions obtained by the iriplicit 
and ID«.C methods. V/o took h = .025 and k = ,01 and/or .1, 
and th j grid was bounded by t = 1. 

Teble 9. 3.1 3rrors by the IDoC ii-^thods with h = ,025 


Probleri 

k 

• 

implicit 

IDeC(l) 

IDeC( 2) 

2 

,01 

7.10-2 

5.7E-2 

5. 43-2 


.1 

1. 2E-0 

l.lB-0 

9.1E-1 

3 






.01 

2.9S-1 

1. 6E-1 

NO 

4 

.01 

H 

\ 

O 

1. 4B-0 

l,0E-0 

9.4 Elli-ptic Partial Differ sntial Equation! 



nox'f consider the Poisson’s equation 
(9.4.1) u^^ + Uyy = g(x,y) 0< x,y <1 

with u(x,y) = f(x, y) on th'^ boundary of the unit squar®. 
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Iri th3 notation of Section 9.2 with t = y, T = 1 
k = h, th'^ five-point finite diff^^renco schvemo for 
( 9. 4. l) is givon by 


‘ 5-'l - 2) ^.1, J 3 + “l, 3+1 + “x, 0-1 - 0 


= h2g(x^,y^), 1,3 = l(l)N-l. 

Th"* above systori of linear eq.uations together with the 
discr'^to boundary conditions is solved iteratively by the 
SOR in'-thod. The iterations were terrn mated when the diffe- 
rence' in two successive iterates is less than TOJi^ 

(iq iKe. ' not'vys 

The NP for (9.4.1) is 


u + u , 
■XX yy 


= P + P 
XX yy 


0 < x^y 4 1 


with the neaning of P as given in Section 9.2. The finite 
difference e^iuations for this NP are the same as (9.4.2) 
but with the righthand side h2[(p^^)^^^ + (P^y)^^^]. 

The approximate values for the dmivatives P^^ and P^^ at 
( X now found uping pi3cewise int srpolatory 

polynomials of degree m which are defined analogous to 
the x-splines and the t-splin'^'s of Section 9.2, 

We took three problems of th° form (9.4.1) with 
f = u end g as given below. Tbe maximum absolute errors 
in th'ir solutions obtained by the five-point scheme and 
th^ IDcC method using univariate polynomials of degree 8 
ar" tebulat-'^d in Table 9.4.1. 



Problo- 1 
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g(x,y) -2[x{l-x) + y(l-y)] 

--xact -olution; u(x,y) = x( l-x)y( 1-y) . 

Iroblan 2 

g(x,y) = 4 /(x 2 +y 2 +i )2 
-'xnct solution: u(x, y) = In(x2+y2+i). 

I’robl^'i 3 

g(x,y) = 0 

exact eolution; u(x,y) = cos ity. 

T able 9.4.1 llrrors by the ^-DOint and IDeC methods 




,1 ..I. I- I..—.,— 

h = 1/8 

h = 1/6 

1 - 

Problen 

TOL 

Pive- 

point 

IDgC(1) 

IDeC(2) 

PlV'^- 

point 

IDeG( 1) 

IDeC( 2) 

1 

.01 

53-4 

73-5 

53-4 

33-4 

43-5 

53- 4 

2 

.ooi; 

43-4 

13-5 

43- 6 

83-5 

13-6 

43-5 

3 

.01 i 

1 

43-2 

13-3 

43-2 

13-2 

83-4 

23-2 


Th^ results show thet there is a slight improveraent 
in the accuracy of the solutions only for a single appli- 
catior of defect corr'^ction. 






5 Conclusion 


In this chnptor we tave introduced a slight rnodi- 
ric''tion in tho Gonstxuctioja of neighbouring problems us(^d 
in the defect correction procedure, and applied th“ modified 
'DO to SO’- classes of problems in partial differential 
equations. 

In th^ modified IDoG applicable to partial differen- 
tial equations wo need to construct global approxireations 
“'S functions of a single variable in contrast to functions 
involving lore th^n one vane blq^ used by Zadunaisky (1976). 

Computational results indicate that the irprovenont 
brought about by the modified IDcC is sometimes substantial, 
and that the first iterate of this IDeC is usually more 
eccurat^ than the other itere^tes. 



GHArTBR 10 


CONCLUSIOMS 


The pruiiary oboective of the work presented in this 
thesis has been to a,pply the method of IDeC in new ways 
and ia new types of operator equations. Consequently, we 
have applied ILeC to integral equations; presented a 
variant of IDeC called Successive Updated IDeC; and tested 
the performance of a modified version of IleC applicable 
to partial differential equations (PDBs), Moreover, we 
have applied ILcC on extrapolation (iDeCB), extrapolation 
on ILeC (SBILeC), ILeC on the deferred correction process 
of Fox (iLEeC), and extrapolation on ILDeC (SBILDeC). We 
have also compared our now methods of iterative improvement 
with other methods of iterative improvement. 

lO * 1 Liscussion of Uumerical Results 
10.1.1 A-pnlicability of ILeC 

Fredholm Integral Equations 

Wo have applied ILeC on quadrature methods for second 
kind FIBs. Wo have considered ILeC^, and ILeCg methods 
based respectively on trapezoid and Simpson's rule* The 
following statements have boon found to be true on the 
basis of our theoretical and numerical results 
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TLj'wrr^ 3.2.1 T,nd tables in Sections 2.2.2 and 3.2). 

( 1) Thr’ higher the order n of the basic quadrature 

' tl-'d, the nore rapid the convergonco of the corresponding 

lx ^C. 

(2) Moro spf'Cif icallj/ , the incretise in order uf 
c>aivc’rg'->ncG uf IloC per iteration is equal to n. However, 
tn naxifiu'" atta^inable order of IleO is independent of n 
but depends on tho dogroe n of polynor-ials used. 

(3) for n given ri, the iterates of ILeC^ and IleCg 
converge to tho sane solution but the latter requires 
fearer iterations then the femer. 

(4) Ihe greater the degree n ( 4 ^6 n 12), the 
hipgher the accuracy in the solution produced by IleC. 

ILoG does not achieve Cin order greater than the 
inaxiL.'UCi attainable stated by the theory, except for a single 
case (Froblem 3 of Section 2.2.2). This is due to the 
fact that the polynomials which extend tne discrete approxi- 
mate solution do not introduce an error due to interpolation 
( since the exact solution of the problem is a polynomial 
of degree one). 

N onlinear Voltorra Integral Equations 

Our co-^putational results for IDeC techniques 11 eO-^, 
IDcC^ and IDeG^^ based on the midpoint rule, the trapezoid 
rule, and Noble's method respectively, show the validity 
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cf the state'^.ents (l) and (4) mentioned above for FIBs. 

Thus IleCjj yields more accurate solution than IDeO^, and 
li/cC^ more accurate solution than IDeCj^. Our IDeC^ method 
turns out to be better than Garey's modified increment 
a.)thudo nnd Wotravali's spline approximation method for 
nonlin^^ar VIFs. 

\!o have also found that for a given stepsixo, the 
11 oC tcchniq.ues generally produce more accurate solution 
than any other discretization method. 

Linear Volterra Integral 'Bg.uations 

The ILoC method based on q.uadrature method that uses 
the trapezoid rule for the solution of linear VIFs produces 
as good results as the IDeC^ method for nonlinear VITIs 
(Tables 2.4.1 and 8.2). 

10.1.2 Annlicabilitv of S^IDeC 

In the case of FI3s, the numerically estimated order 
of the SEILeC method has been found to agree with the 
theoretically predicted order (see equation (3* 3*3) and 
Tables 3.4.1 and 3.4.2). 

10.1.3 Annlicabilitv of IDLeO 

Fredholm Integral Fq.uat;;ons 

The computed orders of convergence for the iteiations 
of the IILeC method based on Gregory's formula for FIFs 
Tvrpp with the theoretical orders (see equation (4.3.19a) 
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'mil T'T.tlco ^.4.4 to 4.4.7), 

t/j mv." also pointod out tha,t IlLeC based on trape- 
ze, id rul.-' WD th correction is easier to mplenent than III eC 
b"'!; “tl jti 3ir<pson*s rule with correction. 

Vult>~rra Integral Equations 

The local approach of ILC has been found to be 
cuporior to the global approach, and IlLeC based on the 
local ILC better than the one based on the global IIC. 
llowcvf'r in the case of nonlinear Vlfls, our results shoiir 
that the ILLeO based on the global IDC yields results 
alr'ost as accurate as the ILIeC based on the local ILC, 
oven thi-iugb the iterates of the global IDC fail to converge 
(see Table 4.4.14b). 

ho marks 

In the case of VIBs, whether linear or nonlinear, 
nunerical results obtained by II eC and IDI eC methods do 
not indicate any asjaaptotic behaviour concerning the order 
of convergence of the methods. 

A comparative study of our methods with certain 
other methods of iterative improvement has been made in 
Chapter 8. 

For linear VIB (2.4.1), we also implemented the ILeC 
which makes use of the following idoa of error estimation 
suggested by Frank ( s<=e Section 6, Frank[1976] ) : discretize 


tho equation 
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X X 

y(x) - k(x,t)y(t)dt = 2f(x) - P„( x) + fkl x, t )r ( t)dt 

Hi w/ n 

0 0 

"tnd takA th'^ results directly as better approximation to 
y(x^). In our experiments, this modification produced the 
eai**'' results as the lueC based on Zadtmaisky ’ s idea, 

10 . 2 Oonclusions 

Our experinents show that 'f«r the construction cf 
p *rturbationo or defects in ILeC, the use of piecewise 
pulynomial interpolation is far more rewarding than that 
of cubic spline interpolation. 

Per solving second kind FIPs, each of the methods 
lueC, SBILeC, II]jeC and SEILLeC is very powerful and 
effective in practice. However, the method IDIeO/SBIlLeC 
based on the trapezoid rule or the method SBILeC based on 
Simpson’s rule is to be chosen for use in general purpose 
routines. 

In the case of general purpose routines for sol\ring 
VIBs, the IliLeC method based on the local approach of II C 
which uses Gregory formula is to he preferred. 

When moderately accurate solutions are desired at 
low costs, our newly introduced methods come in handy. 

In the case of integral equations, standard software 
IS available only for second kind linear FIBs with smooth 
kernels and Green's function kernels (Lelves,' et al[l98l]). 
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But -ur work indiop-tos that the methods such as IleC, S'^ILqC, 
I*^,o0 and SUIjj^-gC will he of use in the development of 
st'indnrd f^-enoral purpose routines for the solution of 
int* /’iral and f'^r other operator equations. 

cni': the methods of iterative improvement, techniques 
based on ll eC are exp'-^cted to he prominent in the present 
tiar* of c”'phasis on mathematical software for large computers. 

10-5 Future Uork 

Finally our work on methods of iterative improvement 
points the wa.y for future research in numerous new areas. 

Some of these are the following. 

(a) I aplementation of ILeC to nonlinear FISs and 
integro-dif f erential equations. 

(h) Applicability of SEILeC to nonlinear FISs and 
ordinary differential equations, 

(c) Implementation of SSIleC on simple numerical 
methods for integro-diff erential equations, with 
the use of the results of extrapolation methods 
(e.g. Chang[l982]). 

(d) Application of ILLeC and SSIIL'oC to ordinary 
and partial differential equations. 

(e) Implementation of IlleC based onlereyra's version 
of deferred correction. 
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(f) i^pplc^ition of SUILeC with the use uf interpo- 
lating polynomials of varying degrees to initial 
value prohlens of ordinary differential equations, 

luring the final stages of our work, we came across 
a ru-w type of iterative improvement technique of BOhmer [I98l] 
calii'd hy him the Eiscreto Newton Method (DNM). This 
peth'id is 'btained as a combination of a Newton and a dis- 
cretization Method. It produces approximations of orders 
2n, 3n, . . . , where n is the order of accuracy of the initial 
numerical solution, BOhmer also presents a general theory 
for INMs, ILeCs and IDCs applicable to any discretization 
method with an casymptotic expansion. DNM is superior to 
and cheaper than IDeC and IDG (Pereyra's version) for 
problems requiring solutions of nonlinear systems of 
equations. 

In future we would like to compare the efficiency of 
our methods with BNMs taking into consideration parameters 
such as GPU time, storage and number of function evaluations. 
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